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ABSTRACT 

Aims. We present a study of protoplanetary disks in spatially resolved low-mass binary stars in the well-known Orion Nebula Cluster 
(ONC) to assess the impact of binarity on the properties of circumstellar disks. This is currently the largest such study in a clustered 
high-stellar-density star-forming environment, as opposed to previous studies, which have mostly been focussed on the young, low- 
, stellar-density Taurus association. We particularly aim to determine the presence of magnetospheric accretion and dust disks for each 

D . binary component, and measure the overall disk frequency. 

Methods. We carried out spatially resolved adaptive-optics-assisted observations to acquire near-infrared photometry and spec- 
troscopy of 26 binaries in the ONC, and determine stellar parameters such as effective temperatures, spectral types, luminosities, 
and masses, as well as accretion properties and near-infrared excesses for the individual binary components. On the basis of our 
medium resolution A"-band spectroscopy, we infer the presence of magnetospheric accretion around each binary component by mea- 
suring the strength of the Brackett-y emission. The accretion disk frequency among the ONC binaries is then estimated from Bayesian 
statistics. The observed disk signatures, measured accretion luminosities, and mass accretion rates are investigated with respect to the 
binary separation, mass ratios, and distance to the center of the ONC. 

Results. A fraction of 40*J°% of the binary components in the sample can be inferred to be TTauri stars possessing an accretion 
disk. This is only marginally smaller than the disk fraction of single stars of ~50% in the ONC. We find that disks in wide binaries of 
>200 AU separation are consistent with random pairing, while the evolution of circumprimary and circumsecondary disks is observed 
to be synchronized in close binaries (separations <200AU). Circumbinary disks appear to be unsuitable to explain this difference. 
Furthermore, we identify several mixed pairs of accreting and non-accreting components, suggesting that these systems are common 



00 



and that there is no preference for either the more or less massive component to evolve faster. The derived accretion luminosities 
and mass accretion rates of the ONC binary components are of similar magnitude as those for both ONC single stars and binaries in 
the Taurus star-forming region. The paper concludes with a discussion of the (presumably weak) connection between the presence of 
£f>l ■ inner accretion disks in young binary systems and the existence of planets in stellar multiples. 

Key words. Stars: late-type - Stars: formation - circumstellar matter - binaries: visual 

1- Introduction for the formation of planetary systems. If disks in binaries 
* \ evolve differently from those in single stars, then the proper- 
One of the most important and explicit indicators of stellar youth ties of the population of plane ts found in binaries can reflect 
^ ■ are circumstellar disks, which seem to be a universal feature those di ff er ences. Today, more than 40 planets are known to 
of the star-forming process. The existence of a stellar compan- orbit one of the compone nts of a binary system, with binary 
ion close to a disk-bearing star, however, has a complex im- separations of mostly >30 AU dMugrauer & Neuhausen f2009: 
pact on the evolution of a disk: it will be truncated to about Eggenberg er & Udrvl2010l) . Interestingly, systems in which both 
a third of the binary separation, may become eccentric and componen ts orbit their own planet are disproportionally rarely 
warped, and the material in the disk may be heated and dy- obserV ed - to the knowledge of the authors no system has been 
■ namically stirred ( Artvmowicz & LubowU994j | Kley et al. | 2008t published so far. Although not free from observational selection 
| Fragner&Nelson| | 2010| ; | Nelson| | 2000|) . Truncation, in particu- effectSi this might indicate that there are differences in the evo- 
lar, can explain observations of young low-mass binary compo- lutions of the initial circums tellar material around the individual 
nents that have a systematically lower dis k frequency than single components of a binary system 
stars with otherwise identical propert ies (iBouwman et al.l l2006: 

iMonin et af I l2007t ICieza et all [2009). Material that is usually Another open issue is universality throughout the Galaxy, 

transported inwards from the outer parts of the disk is missing, Disk evolution is not determined entirely by the interaction 

thus cannot replenish the inner disk and extend their disk life- with their host stars or stellar binary companions. The irradi- 

times, in a way assumed for single star disks. ation by nearby hot sour ces and the gravitation al environment 

Primordial circumstellar disks are not only indicators of play an important role (IScallv & Clarkd feOOll) . While previ- 

the formation of the star itself, but also contain the material ous studies of the signatures of disks in binaries were mostly 

focused on nearb y, loose, stellar associations such as Taurus 

* Based on observations collected at the European Organisation for and Ophiuc hus (IWhite & Ghedl200ltlHartigan & Kenvonl2003: 

Astronomical Research in the Southern Hemisphere, Chile. ESO Data Prat o et alj l2003l iDuchene etalj l2004t iPatience et alj [2008; 
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ID: 074.C-0575. ICieza et aUl2009h . there is evidence that most field stars are in 
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stead born in d ense high-mass clus ters such as the Orion Nebula 

Cluster (ONC: ILada & Ladall2003h 

As t he nearest (414± 7pc; iMenten et ail 120071) young 
(~1 Myr; Hillen b"randl[l997l) high-mass star forming cluster, the 
ONC has been targeted by many studies owing to its abun- 
dant c i rcumstellar d i sk (e.g . IHillenbrand et alJ [T998; Lada et al 



2000; DaRioetal 
1994t iPadgett et al 



20101) and binary content (Prosser et al 
119971: lPetretal.1 1 19981: ISimon etalJll999t 



low spectral resolution to detect spectroscopic pairs. In the fol- 
lowing, we treat these binary components as single stars. 

Table Q] lists all targets including the dates of their ober- 
vation, the instrument used, and the observational mode. 
Furthermore, we provide the projected distance to the massive 
and bright 6 l Ori C system, which is at the center of the ONC. 



KohleretaL|2006; Reipurth et alJ2007l) . although a detailed sur- 2 - 2 - NACO/VLT 



vey of disks around binary components has yet to be performed. 
In this paper, we present the first spatially resolved near-infrared 
spectroscopic observations of a large sample of binaries in a 
clustered star-forming environment - the Orion Nebula Cluster - 
to investigate its circumstellar disk content and the binary com- 
ponent properties. 



2. Observations and data reduction 

2.1. Sample definition 

We ob serve d 20 visual b inaries in the ONC from the binary cen- 
sus of|Pe3 (I1998I) and iKohler etail d2006l) . Data for six addi- 
tional binaries were taken from observations by ICorreia et al.l 
(I20121 in prep.). All targets are listed in Table Q] The projected 
separations range from 0.25 to 1.1 arcsec, which corresponds to 
roughly 100^-00 AU at the distance of the ONC. Magnitude dif- 
ferences of the binary components range from 0. 1 to ~3 mag in 
H and /f s -band. 

All of our targets are members or very likely members 
of the ONC, w hich were mostly identified based on their 
proper motion (IHillenbrand! 1 1 9971 and references therein). 
Targets without proper motion measurements ([HC2000]73, 
TCC 15, TCC55) or those with low proper-motion member- 
ship probability (JW235, JW566, JW876) were confirmed to 
be young stars, thus likely members b ased on their X-ray 
activity (IHillenbrand & Carpenter! I2000I: iReipurth et ail 120071; 
iGetman et al.ll2.Q05h . No further information is available for two 
targets, [AD95] 1468 and [AD95] 2380, but our spectroscopy 
shows late spectral types at moderate luminosity and extinc- 
tion for these and all other targets, ruling out foreground and 
background stars. Common proper motion with the ONC and 
signs of youth combined with small angular separation of the 
components render it likely that all binaries are gravitationally 
bound. However, chance alignment of unrelated members of the 
Orion complex cannot be excluded. A possible example might 
be TCC 15, whose secondary shows almost no photospheric fea- 
tures but strong Bry and He I emission, suggesting a highly 
veiled nature of this component. This makes it a candidate mem- 
ber of the Orion BN/KL region, which is in the line of sight 
of the O NC but slightly furt her away (~450pc) and probably 
younger ( Ment enet al.l 120071) . However, the chance of finding 
an unrelated stellar component at the separation measured for 
this binary (1'.'02) is only W=8% given a proj ected density o f 
-0.03 stars/arcsec 2 in the center of the ONC (IPetret all [1 998). 
In the following, we treat TCC 15 as a physical binary. 

We also searched the literature for possible spectroscopic 
pairs among our binary component s and found that among the 1 1 
binar ies that had been surveyed (Tobin et al. 2009; Fures z et"aT] 
2008), two are spectroscopic binary candidates (see Table []]). 
Both had already been excluded from our statistical analysis for 
other reasons (see discussion in Sect. I3.5.2l for JW 260; JW 974 
was not observed with spectroscopy). Spectroscopic binarity 
cannot be excluded for the other targets, since no published spec- 
troscopic binary surveys exist and our observations are of too 



Twenty of our 26 targets were observed with the NAOS- 
CON ICA instrument (NACO: iLenzen et alOool iRousset et al.l 
2003) at UT4/VLT in the time from November 2004 to February 
2005. Of those, 16 were observed in spectroscopic mode and 1 1 
in CONICA imaging mode (see TableQ]). All observations were 
executed using the adaptive optics system NAOS where the tar- 
gets themselves could be used as natural guide stars. Depending 
on the brightness of a target, either an infrared or visual wave- 
front sensor was used. 



2.2.1. NACO Imaging 

We employed the S13 camera of NACO with a pixel scale of 
13.26 mas/pixel and 13'.'6xl3'.'6 field of view. Each imaged tar- 
get binary was observed in J and H filters and three of these ad- 
ditionally in K s . Observations of the same target in different fil- 
ters were obtained consecutively to minimize the effects of vari- 
ability. The FWHMs of the observations are typically 07075 in /, 
07065 in H, and 07068 in A^-band. All observations were made 
in a two-offset dither pattern with a 5" pointing offset to allow 
for sky subtraction with total integration times varying between 
60 s and 360 s per target and filter, depending on the brightness 
of the target components. 

Imaging data were reduced with custom IRAF and IDL 
routines according to the following procedure. A low ampli- 
tude (< 10ADU), roughly sinusodial horizontal noise pattern- 
probably 50 Hz pick-up noise at read-out-was removed from all 
raw images by subtracting a row-median (omitting the central re- 
gion containing the signal) from each column. We cropped each 
image to a 400x400 pixel (= 573 x 573) subregion of the original 
1 024x1 024 pix preserving all flux from each target multiple. For 
each target and filter there were at least two dithered images that 
undergo the same corrections, which were then used for sky sub- 
traction. Bad pixel correction and flat fielding were applied, the 
latter using lamp flats taken with the NACO internal calibration 
unit. All reduced images per object and filter were then aligned 
and averaged. The fully reduced images in //-band are presented 
in Fig.Q] 



2.2.2. NACO Spectroscopy 

NACO was used in grism-spectroscopy setup using the S27 cam- 
era and an 86 mas wide slit. The wavelength coverage is 2.02- 
2.53 fim with a spectral resolution of /?~1400 and a resolution of 
27 mas/pix in the spatial direction. The slit was aligned with the 
binary separation vector to simultaneously obtain spectra of both 
stellar components. All binaries were observed in an ABBA nod- 
ding pattern with a 12" nod throw, performing several nod cy- 
cles for the faintest targets. Total exposure times per target were 
120-2600 s. Spectroscopic standards of B spectral type were ob- 
served with the same camera setup, close in time and at similar 
airmass to enable us to remove telluric features. These standard 
spectra were reduced and extracted in the same way as the target 
exposures. 
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Table 1. Targets and observations 



dist. to Imaging or Imaging Date 



Name" 


#'OriC ['] 


Obs. with 


Spectroscopy 


Filters* 


SB C 


(UT) 


[AD95] 1468 


7.18 


NACO 


I&S 


JH 


no 


Feb 09, 2005 


[AD95] 2380 


6.92 


NACO 


I&S 


JH 


no 


Jan 06, 2005 


JW235 


6.85 


NACO 


I&S 


JH 


no 


Dec 19, 2004 


JW260 


9.00 


NACO 


I&S 


JH 


SBl 


Dec 19, 2004 


JW519 


0.80 


NACO 


I 


JHK S 




Feb 07 & 09, 2005 


JW553 


0.52 


NACO 


S 






Dec 07, 2004 


JW566 


7.17 


NACO 


I&S 


JH 




Dec 08, 2004 & Feb 17, 2005 


JW598 


0.71 


NACO 


s 






Dec 08, 2004 


JW648 


1.00 


NACO 


I&S 


JHK S 




Feb 07 & 17, 2005 


JW681 


1.25 


NACO 


s 




no 


Jan 06, 2005 


JW687 


2.06 


NACO 


I&S 


JH 




Feb 09, 2005 & Dec 07, 2004 


JW765 


14.09 


NACO 


I 


JH 


no 


Feb 07, 2005 


JW876 


14.45 


NACO 


s 




no 


Jan 08, 2005 


JW959 


7.98 


NACO 


s 




no 


Jan 08, 2005 


JW974 


14.96 


NACO 


1 


JH 


SB2? 


Nov 19, 2004 


[HC2000] 73 


2.30 


NACO 


S 






Jan 05, 2005 


TCC15 


0.50 


NACO 


s 






Dec 29, 2004 


TCC 52 


0.46 


NACO 


s 






Dec 08, 2004 


TCC55 


0.49 


NACO 


S 






Jan 10 & 11,2005 


TCC 97 


0.42 


NACO 


I 


JHK S 




Feb 07, 2005 


JW63 rf 


8.95 


GEMINI 


I&S 


JHKL' 




Feb 16 & 19, 2008 


JW m d 


6.10 


GEMINI 


I&S 


JHKL' 




Feb 17 & 19, 2008 


jw m d 


4.67 


GEMINI 


I&S 


JHKL' 


no 


Feb 24 & Mar 7, 2008 


JW391 <i 


2.81 


GEMINI 


I&S 


JHKL' 




Feb 19, 2008 


JW709 rf 


3.53 


GEMINI 


I&S 


JHKL' 




Feb 18 & 20, 2008 


JW 867 rf 


5.79 


GEMINI 


I&S 


JHKL' 


no 


Feb 23, 2008 



Notes. (a) Identifiers are suitable for us e with the simbad database (htt p : //simbad . harvard . edu/ simbad/ I. References: AD95: lAfi & Depovl 
119951 : 7W: lJones & Walkerll988l : rCC: lMcCaugh rean & Stauffer l 19941 : HC2000: Hillenbrand & Carpenter 2000. {b) Photometry for this study was 
obtained in these filters and has been listed in Table [3] If fewer filters than the complete set of JHK S has been ind i cated here, photometry ha s 
been ad ded to Table [3] from literature sources when available. (r) Spectroscopic binary status. References: iTobin et alj J2009): Furesz et al. ( 2008). 
From Correia et al. ( 2012', in prep.). 



Spectroscopic exposures were reduced with IDL and IRAF 
routines including fiat fielding, sky subtraction, and bad pixel 
removal. Lamp flats taken with an internal flat screen were fit- 
ted along the dispersion axis and divided. The extraction and 
wavelength calibration of the reduced target and standard ex- 
posures were performed using the apextract and dispcor pack- 
ages in IRAF. The two-dimensional spectrum was traced with a 
fourth order polynomial and extracted through averaging over 
a 10 pixel-wide (~0727) aperture centered on the ~4 pixel- 
wide (FWHM) trace of each target component. The subsequent 
wavelength calibration uses exposures of an argon arc lamp ex- 
tracted in the same traces. All nodding exposures of the same 
target component were then aligned and averaged. The final ex- 
tracted and calibrated spectra span a wavelength range of 20320- 
25440 A with a resolution of 5 A/pixel. 

To allow for accurate telluric-line removal, intrinsic spectral 
features of the standard stars had to be removed. These were 
the Brackett-y (Bry, 21665 A) and the He I (21 126 A) absorption 
lines, where the latter was only observed for spectral types B0 
and B 1 . The removal of in particular the Bry line is crucial be- 
cause we aim to measure the equivalent widths of Bry emission 
in the target spectra. Therefore, the Bry line was carefully mod- 
eled according to the following scheme. Two telluric lines blend 
into the blue and red wings of Bry at the spectral resolution of 
our observations. We modeled these tellurics by averaging the 
fluxes in the Bry region (before telluric removal) of those target 
component spectra that showed neither absorption nor emission 



in Bry. The telluric standard flux was divided by the thus gen- 
erated local telluric model, allowing us to fit the remaining Bry 
absorption line with a Moffat-profile and remove its signature. 
The cleaned standard spectra were then divided by a black body 
curve of temperature r e ff according to their spectral type dCoxl 
2000) to obtain a pure telluric spectrum, convolved with the in- 
strumental response of NACO. 

We observed a wavelength-dependent mismatch of the wave- 
length calibrations output from IRAF/dispcor between the target 
stars and the corresponding telluric standards of up to about one 
pixel. To guarantee a good match of the positions of the telluric 
features, we used the tellurics themselves to fine-tune the wave- 
length calibration of the standard spectra. A customized IDL rou- 
tine computes the local wavelength difference by means of cross 
correlation and corrects the mismatch accordingly. Since the tel- 
lurics are the strongest features in all our spectra, this method 
results in a good match between the telluric features of the target 
and reference, and telluric features could be removed reliably. 

Flux uncertainties in the derived spectra as needed for the;^ 2 
minimization method described in jQ.3l were estimated from the 
fully reduced and extracted target spectra by performing local 
computation of standard deviations. The final reduced and ex- 
tracted component spectra contain a noisy (signal-to-noise ratio 
(S/N) smaller than 20/pix) region redwards of ~ 25000 A caused 
by the low atmospheric transmission at these wavelengths. The 
region with A > 25120 A was excluded from any further evalu- 
ation. The set of reduced and extracted spectra consists of one 
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[AD95]1468 
• 

1 


[AD95]2380 


JW235 

♦ 

• 


JW260 

• 

• 


JW519 

# 

• 


JW566 

* 

• 


JW648 

• 

* 


JW687 

4 

• 


JW765 

- 

• 

N 


JW974 

* 

1 


TCC97 





Fig. 1. H-band images of all targets observed with NACO imaging. The intensity scale is linear and adapted to best depict both 
components. 



spectrum for each target component in the spectral range of 
20320 A - 25120 A with a resolution of 5 A/pixel. All final re- 
duced spectra are displayed in Fig. [2] The positions of the most 
prominent absorption and emission features are overplotted and 
listed in Table [2] 

2.3. NIFS-NIRI/Gemini North imaging and spectroscopy 

The observations with NIRI photo metry and NIFS spe ctroscopy 
at Gemini North are described in lCorreia et alj d2012l in prep.). 
We used the reduced and extracted but non-telluric corrected 
spectra of Correia et al. (2011) and performed the telluric cor- 
rection using the same method as for our NACO observations. 
Furthermore, since the NIFS spectra are of higher spectral reso- 
lution (R~5000) than both our NACO observations and the tem- 
plate spectra, we smoothed the NIFS spectra with a Gaussian 
kernel to a resolution of R~1400. These steps were included to 
guarantee that a coherent evaluation of all data is possible. The 
six target spectra observe d with NIFS that wer e included in this 
study will be presented in lCorreia et al ] (I2012L in prep.). 

3. Results 

3. 1 . Photometry and astrometry 

Relative aperture photometry for all targets observed with 
NACO was obtained by applying the PHOT task in the IRAF 
DAOPHOT package to each of the reduced binary star images. 
The aperture radius was varied from 2 to 20 pixels to find a possi- 
ble convergence of the magnitude difference of primary and sec- 



ondary. The differential photometry of most binaries converged 
for aperture sizes of 3 to 6 pixels allowing determination of their 
magnitude differences with uncertainties of Amag < 0.03. For 
the binaries that did not converge but followed a monotonic de- 
crease in their magnitude difference with aperture size, we as- 
signed a value of Amag by averaging the results for apertures 
of sizes between 3 and 6 pixels. The uncertainty was estimated 
according to the slope of each individual curve. 

The DAOFIND task in DAOPHOT returned sufficiently ac- 
curate astrometry for all targets. The pixel data were transformed 
into physical angles and separations using the pixel scale of 
0.01 3260 arcsec/pixel and the rotation offset of 0°. The result- 
ing relative photometry and astrometry are listed in Table [3] 

The procedure was successful for all binaries observed with 
NACO, although we found peculiarities for two targets. The 
image of the brighter component of JW519 has an elongated 
shape (see Fig.[TJ that is not seen in the companion point spread 
function (PSF). We therefore conclude that we have found evi- 
dence of an unresolved third component in the JW519 system. 
Unfortunately, the close separation of the system does not allow 
us to determine any of the separate parameters of the individual 
components. The photometry of JW519 in Tables [3] and |5J as 
well as all further evaluation, therefore treats this likely triple as 
a binary. The primary of TCC 97 is s urroun ded by a proplyd that 
was first identified by |Q' Dell etall d!993l) . Since this feature is 
detected in our images, our aperture photometry averages only 
the smallest useful apertures of 3-5 pixel radius (instead of 3-6) 
to exclude as much of the disk flux as possible. 

A special photometry routine was applied to the target bina- 
ries JW 553, [HC2000] 73, and TCC 55, which used photometry 
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I I I I I I 
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C0(6V=2) 




JW260 



CO(i5i/ = 2) 



I I I 1 I 1 

JW566 



CO(i5u = 2) 




I I I 1 I 1 

JW648 - 



CO(i5u=2) 





21000 22000 23000 24000 21000 22000 23000 24000 
Wavelength A (A) Wavelength A (A) 

Fig. 2. Spectra of the primary (top spectrum in each panel) and secondary compone nt (bottom spectrum) of all targets observed 
with NACO spectroscopy (Gemini/NIFS spectra are displayed in lCorreia et al.l (I20121 in prep.)). Primary spectra are normalized at 
2.2 jjm, and the secondaries are arbitrarily offset. The position of the most prominent lines in Table|2]are indicated. 
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Fig. 2. Continued 
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Table 3. Relative photometry and astrometry of the observed binaries 1 





A / 

11. 1 


Hit 


A k' 




sep 


PA ' J 




Name 


[mag] 


[mag] 


[mag] 


[mag] 


L J 


r° l 
L J 


Kei. 


[AD95] 1468 


0.75+0.10 


0.66 ±0.02 


0.10 ±0.05 




1.08 


76.9 


T,l 


[AD95] 2380 


>2.5 C 


2.62 ± 0.03 


2.99 ±0.15 




0.59 


77.6 


T,l 


JW235 


0.46 ±0.02 


0.10 ±0.02 


0.47 ±0.15 




0.35 


163.6 


T,l 


JW260 


0.53 ±0.02 


0.40 ± 0.02 


0.17 ±0.05 




0.35 


292.2 


T.l 


JW519 f 


2.8 ±0.2 


2.8 ±0.2 


2.6±0.1 




0.36 


204.3 


T 


JW553 


2.1 ±0.3 


3.2 ±0.3 


3.19±0.10 




0.384±0.004 


248.1±0.3 


2 


JW566 


0.20 ±0.03 


0.39 ±0.02 


0.78 ±0.02 




0.86 


33.8 


5 


JW598 








3.19 


0.9 




4 


JW648 


1.12 ±0.03 


1.16 ±0.03 


1.23 ±0.02 




0.68 


278.7 


T 


JW681 








1.53 


1.09 


214 


3,4 


JW687 


1.63 ±0.05 


1.15±0.07 d 


0.53 ±0.07 




0.47 


232.6 


T 


JW765 


0.08 ±0.02 


0.12 ±0.01 


0.06 ±0.10 




0.33 


16.5 


T.l 


JW876 






0.50 ±0.05 




0.49 




1 


JW959 






0.07 ±0.03 




0.34 




1 


JW974 


1.16 ±0.04 


1.26 ±0.04 


1.41 ±0.05 




0.32 


128.7 


T.l 


[HC2000] 73 


0.48 ±0.09 


1.13 ±0.06 




1.52 ±0.07 


0.71 


266.2±0.5 


4 


TCC15 


2.5 ±0.3 


3.3 ±0.3 


3.65 ±0.1 




1.022±0.004 


288.2±0.3 


2 


TCC 52 


1.61 ±0.03 


1.56 ±0.03 


1.64 ±0.02 




0.52 


39.24 


5 


TCC55 


1.0±0.U 


2.2 ±0.3 


1.26 ±0.10 




0.256±0.004 


153.1+0.3 


2 


TCC 97 


1.8 + 0.15 


1.8 ± 0.15 


1.4 ±0.2 




0.88 


98.5 


T 



Notes. m Repr oduced here are onl y the targets observed with NACO. The photometry of the six additional targets observed with Gemini/NIRI 
can be found in ICorreia et al1 l l2012l . in prep.). (n> Uncertainty in the separation: Asep = C.'Ol unless otherwise noted. ib) Position angle uncertainty: 
APA = 0.5° unless otherwise noted. (c) The companion to [AD95] 2380A is detected with less than 3<r significance in i-band. The number quoted 
here is a lower limit. (< " This number is an average of two independent measurements that differ by 0.16 magnitudes. w The elongated shape of 
JW519's primary suggests the primary to be binary itself. The photometry for JW519, howev er, is based on the assumption of a single central 
object since separate components cannot be identified. (rt Reanalysis of MAD data published in B ouv et aTJ (12008). 

References. (T) This paper; (1) iKohler etail 120061: (2) iBouv et ail 120081 ; (3) iReipurth et alj 120071 ; (4) IPehl [19981; (5) ESO archival data 
074.C-0637(A). Reduction and photometry are described in Sects. [2~2~T1 and [3~T1 



of several reference stars in the same exposure. In particular for 
JW 553 and TCC 55, this procedure leads to more accurate pho- 
tometry in the environment of strong nebulosity and high stellar 
density close to the cluster center. The photometry of JW553 
and TCC 55 uses fu lly reduced J, H, and A'-band mosaics of the 
Trapezium region (Bouv et al. 2008) observed with the MCAO 
demonstrator MAD (Marchetti et al. 2007). A reference PSF was 
computed from 4-8 stars within ~9" of the target binary. With 
this PSF, we obtained instrumental photometry for all stars in the 
vicinity of the target using the IRAF daophot package. Observed 
apparent magnitudes of the reference stars together with our 
measured relative magnitudes were then used to derive the ap- 
parent magnitudes of the target binary components. Similarly, 
PSF photo metry of [HC2000] 73 was derived using ADONIS 
dPetrll 19981) exposures making use of the PSF and appa rent mag- 
nitude of the nearby star 6> 2 Ori A dMuench et al.ll2002l) . 

Most photometry in this paper was obtained in the NACO 
and 2MASS JHK S filter systems. Since these two filter sets are 
very similar, no transformations had to be applied. Some pho- 
tometry, however, was performed in different filters. We checked 
the compatibility of the 2MASS K s and the ADONIS and MAD 
K filters used for some of the systems and relative magnitudes 
(see Tables [3] and @). Comparing 39 young stars in the ONC 
(with similar properties as the target sample) for which both 
ADONIS K and K s photometry are available dPetrlll998l) . we 
found an average offset of only 0.01 +0.23 mag between K and 
K s . The Gemini/NIRI A'-band photometry of late-type stars is 
compatible with that o f the AVband with K - K s = + 0.05 mag 
dDaemgen et al.l 120071) . We therefore did not apply any correc- 
tions to our A'-band magnitudes. Some of the relations used in 



Sect. 13.21 (i.e. the dwarf, giant, and CTT locuses) were trans- 
forme d to the 2MASS system using the relations in ICarpenterl 
d200lt) . In the following, we thus assume that all photometry is 
compatible with the 2MASS filter system. 

The relative photometric and astrometric results are given in 
Table [3] Together with literature values for the integrated pho- 
tometry in Table |4] component magnitudes were derived, which 
are listed in Table [5] 



3.2. Color-color diagram, extinctions, and color-excess 

Using the component magnitudes of Table [5] we composed 
a (H-K S )-(J-H) color-color diagram (Fig. [3}. We compared 
our d ata with the loci of dwarfs and giants (Bess ell & Bretll 
1988) and the location of classical T Tauri stars (CTTs lo- 



cus; iMeveret all 1 1997b that had both been converted to the 
2MASS photometric system. Extinctions were derived by dered- 
denin g to the CTTs locus along the interstellar reddening vector 
dCohen et al.ll 198 lb and are listed in Table [6] 

Most targets are located in the region accessible from the 
dwarf and CTTs loci by adding extinction. However, three 
groups of targets have rather peculiar locations in the diagram: 
i) two of the targets with among the smallest of H-K s val- 
ues, the secondaries of JW63 and JW 176, have no intersection 
with either the dwarf or giant locus along the dereddening di- 
rection, ii) There are some targets (the secondaries to JW553, 
[HC2000] 73, and TCC 97) that are significantly below the CTTs 
locus. Targets if this location in the col or-color diagram hav e 
been observed before (see e.g. Fig. 20 in lRobberto et alJ l2010). 
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Table 4. Non-resolved photometry of the observed binaries 





y S y S 




K? 


s 


K sys 




Name 


[mag] 


[mag] 


[mag] 


[mag] 


Ref. 


[AD95] 1468 


13.92 + 0.11 


12.48 ±0.18 


1 1.23 ± 


0.08 




3" 


[AD95] 2380 


13.9 ±0.02 


11.39 ±0.03 


9.8S ± 


U.U2 




3" 


JW235 


12.10±0.15 


11.11 ±0.15 


10.50 ± 


0.15 




2 


JW260 


8.19 + 0.15 


7.60 ±0.15 


1 .Li ± 


U. 13 




2 


JW519 


12.07 ±0.01 


11.07 ±0.01 


lU.Ol ± 


AAI 

U.U1 




1 


JW553 


10 56 + 10 


9.49 ± 0.05 






9 05 +0 10 


b 


JW566 


11.49 ±0.04 


9.97 ±0.05 


O Q£. i 

8.80 ± 


U.U3 




3" 


JW598 


10.87 ±0.04 


9.57 ±0.05 


o ni i 
8.93 ± 


U. 11 




2 


JW648 


10.97 ±0.03 


9.91 ±0.02 


9.28 ± 


0.01 




1 


JW681 


12.78 ±0.01 


11.82±0.01 


1 1 .07 ± 


0.11 




1 


JW687 


11.94 ±0.01 


10.50 ±0.01 


9.70 ± 


0.03 




1 


JW765 


11.76±0.15 


11.04 ±0.15 


10.81 ± 


0.15 




2 


JW876 


9.31 ±0.02 


8.45 ±0.03 


8.03 ± 


0.02 




3" 


JW959 


9.36 ±0.02 


8.84 ±0.03 


8.63 ± 


0.02 




3" 


JW974 


12.42 ±0.02 


11.78 ±0.03 


11.41 ± 


0.02 




3" 


[HC2000] 73 


12.59 ±0.05 


11.72 ±0.03 


10.99 ± 


0.03 




b 


TCC15 


12.96 ±0.03 


11.14 ±0.01 


10.25 ± 


0.01 




1 


TCC52 


8.64 ±0.04 


7.56 ±0.04 


6.72 ± 


0.04 




2 


TCC55 


15.10 + 0.14 


13.27 ±0.14 






11.15 ±0.14 


b 


TCC 97 


13.13 ±0.03 


12.44 ±0.02 


11.77 ± 


0.01 




1 



Notes. <n) If no other reference could be found and if the distance to 1 OriC is larger than 5', 2MASS values were used. (,,) System magnitudes 
were derived from the component magnitudes (Table|5](. 

References. (D lMuench et all 20021 ; (2) ICarpenter et alj2001l ; (3) 2MASS. ICut"ri et alj|200l . 



Table 2. Spectral features identified in the observed spectra 



,1, 

[A] 


Width 
[A] 


Species 


Transition 


20338.0 




H 2 


v= (1-0)5(2) 


20587.0 




He I 


2p'P°-2i 1 S 


21066.6 




Mgl 


r 2,3,4 'S u 3,4,5 


21098.8 




All 


4p 2 P° /2 -5s 2 S [/2 


21169.6 




All 


4p 2 P° /2 -5s 2 S l/2 


21218.0 




H 2 


y = (l-0)S(l) 


21661.2 


56 


HI 


n = 7 - 4 (Bry) 


21785.7 




Si I 




21825.7 




Si I 




21903.4 




Til 


a 5 P 2 -z 5 D° 


22062.4 


116 


Nal 


4p z P° 3n -4s 2 S U 2 


22089.7 


* 


Nal 


4p 2 P° l/2 -A S 2 S l/2 


22614.1 


91 


Cal 




22631.1 


* 


Cal 


4fF°-4d 3 D 2 


22657.3 


* 


Cal 


4/ 3 F°-4rf 3 D, 


22814.1 




Mgl 


4</7), ;j 6/7 V u 


22935.3 


170 


12 CO 


v = (2-0) band head 


23226.9 


170 


12 co 


y = (3-1) band head 


23354.8 




Nal 


4p 2 P° l/2 -Ad 2 D y2 


23385.5 




Nal 


4p 2 P° V2 -4d 2 D 5/2 


23524.6 


170 


12 CO 


v = (4-2) band head 


23829.5 


170 


12 co 


y = (5-3) band head 



Notes. Identified features in our A'-band spectra. Integration widths are 
listed here for all lines that have equivalent widths measured in this 
paper. Lines marked with asterisks blend into the line with the next 
shortest wavelength (the respective previous line in the list); these lines 
were integrated to geth er. The transition informa tion is composed from 
iPrato et all |2003) and Kleinmann ~& Halll dl986h . 




i 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
H-K s (mag) 

Fig. 3. Color-col or diagram of the target binary components. 
The CTTs locus ([Meyer et al.ll 19971) . the dwarf and giant locus 
dBessell & B rett 1988)), and a reddening vector of 10 mag length 
are o verplotted, after conversion to the 2MASS photometric sys- 
tem (ICarpenterl200 lh . Targets to the right of the dash-dotted line 
have an IR color excess. 



Hi) The outlier TCC 55B in the bottom-right of the plot. Most of 
these peculiar locations can be explained by the intrinsic photo- 
metri c variability of young stars of ~0.2mag ( Car penter et al.l 
1200 lh and the fact that photometry was not taken simultane- 
ously. Furthermore, color s are known to depen d on the incli- 
nation of a possible disk (Robita iTie et alj |2006) - a parameter 
that cannot be determined with our data. All peculiar objects 
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Table 5. Individual component apparent magnitudes ' 





J 


H 




Name 


[mag] 


[mag] 


[mag] 



r a r^n^ n 1 a ao 
[AJJVjJ 14oo 


A 

A 


1 A 

14. 


3 A 
JU 


± U. 


1 1 

. 1 1 


i o 
12 


AC 

V3 


+ 


A 

u. 


1 Q 

. 1 5 


1 1 
1 1 . 


yj 


+ 


A 

u. 


AC 
.1)6 




R 


1 5 

1 J . 


1 1 

. 1 1 


m w. 


1 1 


1 3 

1 J 


(SI 




n 


. 1 


1 7 

1 z. 


01 




A 
\t 


08 


r a r^n^n o-ion 
[AlJbDJ ZjoU 


A 

A 










1 1 

1 1 


.45 


+ 


A 

u 


A3 


A 

y. 


AC 

y3 


+ 


A 

u 


AO 

.1)2 




R 










1 4 


1 

. i w 




n 
\t 


04 

AM 


1 7 

1Z. 


04 
y~f 




A 
\t 


1 4 


T\l 7 TIC 

J W Idd 


A 

A 


1 O 

12. 


AC 

.03 


± U. 


1 c 

.13 


1 i 
1 1 


O 1 

.51 


± 


A 


.13 


1 1 

1 1 . 


A 1 

04 


± 


A 

u 


.16 




R 


1 3 

1 J, 


1 1 

. 1 1 




1 J 


1 1 
1 1 


Q 1 
,y i 




A 
\t 


1 c 

. 1 J 


1 1 

1 1 . 


C 1 

. j i 




A 

U. 


1 8 
. 1 o 


JWT O/iA 


A 

A 


5. 


. / 1 


_i_ a 
± U. 


1 c 

.13 


Q 

5 


1 "7 
. 1 / 


+ 


A 
I). 


.13 


n 

1 . 


AA 


+ 


A 

u. 


1 c 

.13 




r 


q 


.Z4 


m w 


1 c 

. 1 J 


8 
o 


S7 




A 
v> 


1 c 

. 1 J 


C 

o. 


07 

M / 




A 
U 


1 c 
. i j 


JWT z i n 


A 

A 


i o 


.13 


i A 

± U. 


.02 


1 1 

1 1 


i c 
.13 


+ 


U 


A") 

,U2 


1 A 


1 A 

. /o 


+ 


A 



A 1 
.01 




R 


1 4 


OS 

. 7 J 


+ o 

it: w. 


1 Q 

. i y 


1 3 


QS 

. jj 




A 
W 


1 


1 1 
i j . 


10 

. J7V7 




A 

U. 


00 

,W 7 


JWT ^^^2 


A 

A 


i a 
IU. 


AA 

.00 


_i_ A 
± U. 


i o 
. 12 


A 

y 


.33 


+ 


A 

u 


.1)0 


A 


. 1 1 


+ 


A 
I) 


.U4 




R 

D 


1 1 


1 8 
. 1 o 


4- 

it: w. 


HQ 
\>y 


1 9 
1 z 


70 




A 


OA 


1 7 

1 z. 


7 1 

,Z 1 




A 

U. 


1 1" 


jwt c/;/; 

J W ODD 


A 

A 


1 o 
12 


.13 


_i_ A 
± U. 


A-1 

.04 


1 A 
1U 


C 1 

.34 


+ 


A 

u 


AC 

,U3 


A 


OA 

,2y 


+ 


A 



A3 

.03 




R 


1 ? 
1 z 


IS 


+ o 

IE W 


04 

,U4 


1 
1 w 


. jj 




A 
v> 


OS 


1 


07 

M / 




A 

U. 


01 


JWT A/1 Q 

J W o4o 


A 

A 


i i 
1 1 


.31) 


_i_ A 
± U. 


A3 

.1)3 


1 A 

1U 


./J 


+ 


ft 

u 


A") 

,U2 


A 


.35 


+ 


A 

L) 


A 1 
.01 




R 


1 7 
1 z. 


47 


4- 

it: w. 


04. 
W4 


1 1 

1 1 


. J 7 




A 


01 

. V7J7 


1 


8 1 

,o 1 




A 

U. 


07 
,WZ 


JWT AQ7 
J W Do / 


A 

A 


1 o 
12. 


1 A 

.10 


_i_ A 
± U. 


A 1 
All 


1 A 
10 


.52 


+ 


A 

u 


AO 

.02 


1 A 

10 


oo 
.22 


+ 


A 



A 1 

.04 




R 




70 
. / y 


+ o 

it: w. 


04 

W4 


1 1 
1 1 


07 




A 
\f 


OS 


1 


oc 
. / j 




A 
w 


OS 


JWT HCQ. 


A 

A 


i o 

12. 


.4/ 


i A 

± U. 


.13 


1 t 
1 1 




+ 


A 

u 


i c 
. 13 


1 1 

1 1 


.30 


+ 


A 



,lo 




R 


1 7 
1 z 


cc 

. j j 


+ o 

IE w 


1 <v 

. 1 J 


1 1 
1 1 


8S 

. 5 J 




A 
v> 


1 C 

. 1 J 


1 1 

1 1 . 


SO 




A 
U 


1 A 


JWT Q"7A 

JWo/O 


A 

A 




















o 
O 


CA 

.30 


+ 


A 



A3 

.03 




R 




















q 

7. 


06 

.WO 




A 
U 


04 


jwt ocn 
J W yjy 


A 

A 




















A 

y. 


i c 
.33 


+ 


A 



AO 

.02 




R 




















q 

7. 


47 




A 

U. 


01 


jwt an a 
J W y /4 


A 

A 


1 o 
12. 


"7 -1 


_i_ A 
± U. 


AO 

.112 


1 O 

12 


no 
.1)0 


+ 


A 

u 


A3 
.1)3 


1 1 

1 1 


.0 / 


+ 


A 
I) 


AO 

,02 




B 


13. 


.90 


±0. 


.04 


13 


.34 


+ 





.04 


13 


.08 


+ 





.04 


[HC2000] 73 


A 


13 


.13 


±0 


.06 


12 


.05 


± 





.03 


11. 


.23 


+ 





.03 




B 


13. 


.61 


±0. 


.07 


13 


.18 


± 





.05 


12. 


.75 


+ 





.06 


TCC15 


A 


13. 


.06 


±0 


.04 


11 


.19 


± 





.02 


10. 


.29 


+ 





.01 




B 


15. 


.56 


±0. 


.27 


14 


.49 


+ 





.29 


13 


.94 


+ 





.10 


TCC52 


A 


8 


.86 


±0. 


.04 


7. 


.79 


+ 





.04 


6 


.94 


+ 





.04 




B 


10. 


.47 


±0. 


.05 


9 


.35 


+ 





.05 


8 


.58 


+ 





,04 


TCC55 


A 


14. 


.45 


±0. 


.08 


12 


.67 


+ 





.17 


10 


.96 


+ 





.03° 




B 


15. 


.45 


±0 


.10 


14 


.84 


± 





.31 


12. 


.22 


+ 





.08" 


TCC 97 


A 


13. 


.32 


±0. 


.04 


12 


.63 


± 





.03 


12 


.03 


+ 





.04 




B 


15. 


.12 


±0 


.13 


14 


.43 


± 





,13 


13 


.43 


+ 





.16 



Notes. The photometry of the six targets observed with Gemini/NIRI 
is listed in ICorreia et al.1 d2012l in prep.). (a) These values are derived 
from Zf-band system magnitudes instead of K s (see Table |4}. 



a substantially larger Ay from the optical data than from our NIR 
measurements. This discrepancy for this target, however, could 
be explained by its large veiling value (optical veiling can reduce 
the measured V-I color and thus the derived extinction). 



3.3. Spectral types and veiling 

Owing to the young age of the targets in our sample, the stars are 
still contracting, hence their surface gravities (logg) are lower 
than those of main sequence stars. However, to find appropriate 
templates for spectral classification and the estimation of both 
the visual extinction Ay and near-infrared continuum excess in 
/T-band (veiling, r^), one needs to find templates with phys- 
ical conditions as close as possible to the pre-main sequence 
stars in this sample. Since no comprehensive catalog of pre- 
main sequence spectra at a spectral resolution of R = 1400 or 
higher exists, we compared our target spectra to thos e of dwarfs 
and gi ants using a method similar to that shown in iPrato et all 
(2003|). We measured the equivalent widths of the T e g and 
log ^-sensitive photospheric features Nal, Cal, 12 CO(2-0), and 
12 CO(4-2) (wavelengths see Table |2| for all target component 
spectra as well as dwarf and giant spectra from th e IRTF spec- 
tral libary (Ra vner et aHl2009{ICushing et al.ll2005l) in a spectral 
range between F2 and M9. Fig.[4]demonstrates that dwarf spec- 
tra are more suitable templates for our sample than giant stars in 
the same spectral range. We note that veiling reduces the equiva- 
lent widths of spectral lines according to w™ easured = Wa/(1 +r K ) 
(see also Sect. 13.5.11 ) and points in Fig. [4] would move towards 
the origin if the veiling is high. It would therefore be necessary to 
only compare targets with small r^, which is only possible after 
the veiling was determined using the templates. However, we see 
that not only the distribution of < 0.2 targets but all targets 
are well-congruent with the dwarf locus, indicating that dwarfs 
are suitable for the determination of the stellar parameters of our 
pre-main sequence stars. We also considered intermediate solu- 
tions between dwarfs and giants as possible templates. Despite 
the possibly closer match in logg of pre-main sequence stars, 
no improvement in the resulting match with our target equiva- 
lent widths could be seen. Thus, in order to minimize additional 
noise sources, we used dwarf templates for the evaluation that 
we now describe. 

Spectral types, spectroscopic extinctions, and veiling were 
simultaneously determined by a x 1 minimization method that 
modifies the template spectra according to 



were assigned an ext i nction of 0. Using expected dwarf colors 
from lBessell & Breta d 19881) . we then derived the color excesses 
E j-h — (J - H)obs - (J - H)o and Eh-k^ for all objects (Tabled 
We compared our exti nction measurements with optical data 
from iProsser et alj d 1994ft . Their sample contains six spatially 
resolved binaries that are also part of our sample (JW553, 
JW598, JW648, JW681, JW687, TCC 153) ■ Using the ob- 
served V-I colors and an estimate of (V-/)o derived from our 
measured effe ctive temperatures a nd the lMyr stellar evolution- 
ary model o f iBaraffe et ail (f 1998), we calculated extinctions as 
A v * \(V - 1) - ( V - 7) n ] /0.4 dBessell & Brettll 1988b . The results 
are consistent with our near-infrared extinction measurements 
usually to within ~ 1 .0 mag, except for JW 687 A, where we find 

1 IProsser et al.1 i 19941) list the brighter component in V as the primary 
of binary. We show, however, that it is of later spectral type than its 
companion. Accordingly, o ur designation swaps both components with 
respect to the lProsser et al.1 paper. 



(cf. IPrato et ai1l200l) with t a = (0.522/A) 1 - 6 A s v pec , where F phsph 
is the photospheric flux of the templates and Ic is the K-bmd 
excess in units of F p h sp h(2.2//m)/c, i.e. the excess over the pho- 
tospheric flux at 2.2yi/m normalized with a constant c. We intro- 
duce the extinction variable ec , since it is typically not identi- 
cal to the photometric extinction Ay. The excess k is assumed to 
not vary strongly with wavelength within the limits of /T-band, 
i.e. it is k(A) - const. Although some of the targets show evi- 
dence of a slightly stronger excess towards the red edge of the 
/T-band, the theoretically slightly poorer fit of line depth at the 
long wavelength end does not have a strong impact on the result- 
ing k since it is determined from the best fit to the entire wave- 
length range of our spectra. For each template spectral type, the 
three variables A^ ec , k, and c were modified with 120-160 steps 
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i i i i 




W x (CO[2-0] + CO[4-2]) 

Fig. 4. Top: The equivalent widths of Nal+Cal versus CO[2- 
0]+CO[4-2] as an indicator of T e g and logg for dwarf (open 
circles) and giant (filled circles) template spectra from the IRTF 
spectral library. Spectral types are color-coded according to the 
legend in the upper left. The red and blue lines guide the eye to 
the two relations. Bottom: The same plot including the derived 
dwarf and giant relations (red and blue curve) as in the top panel, 
overplot with the ONC target components. Open symbols show 
targetsted with high veiling r# > 0.2, filled symbols rx < 0.2. 
Two targets, TCC 15B and [AD95] 2380B, are off bounds at high 
^(Nal+Cal). 

each within a reasonable range of values and we found a mini- 
mum value of 

X n-dof-1 4-* AF 2 

where n is the number of pixels in the spectrum, dof = 3 the 
degrees of freedom, and Fj, AFj, and F* the flux in the 2-th pixel 
of the target spectrum, its measurement error, and the modified 
template spectrum, respectively. The minimum of the x 2 distri- 
bution was compared for different spectral types and the nine 
best-fitting solutions (spectral type along with the correspond- 
ing optimal combination of Ay ec , k, c) when examined by eye. 
Spectral types were selected by comparing the overall shape of 
the continuum and the strength of several photospheric absorp- 
tion lines of the modified (eq.Q]) best-fit templates. Uncertainties 
in the spectral type were estimated from the range of spectra that 
could possibly fit the data. This resulted in a typical uncertainty 
of one or two subclasses. The uncertainty in k was determined 




0.0 r i i i , . .1 

21000 22000 23000 24000 
Wavelength A [A] 



Fig. 5. The result of the spectral template fitting for the primary 
and secondary component of JW 876. The black curves show the 
spectrum of the primary and secondary component respectively, 
offset for clearer visibility. The red curves are the corresponding 
best-fit models from the IRTF spectral library, modified accord- 
ing to eq. (Q]i with extinction and veiling values from Table [6] 

from the e rror ellipse in a three-parameter ^-minimization at 
xLn + 3 5 <tWalJlLl996|). 

An example of the best-fit result for one of our target binaries 
is shown in Fig. [3] From the excess flux k, the /f-band continuum 
excess is calculated to be 



■F p hsph(2.2/mi)/c 

The results from the spectral fitting are summarized in Table [6] 
Some spectra do not show any photospheric features and no es- 
timation of the spectral type was possible. These target com- 
ponents are marked with ellipses in the spectral type-dependent 
columns of Table [6] 

Despite the wide range of extinction values derived from the 
photometric and spectroscopic determination, we did not force 
the x 1 minimization to match the photometric Ay values from 
Sect. 13.21 for the following reason. In the x 2 fitting routine, the 
determination of the other parameters (spectral type, veiling, 
normalizing constant c) does not depend on the value chosen 
for Ay as long as it has reasonable values. While spectral types 
are mainly determined by line ratios, veiling is sensitive to line 
depths; both signatures are not strongly influenced by extinc- 
tion. However, forcing Ay to a particular value typically pro- 
duces poor fits and hence Ay is kept as a free parameter. 

We only use the photometric extinctions listed in Table [6] for 
further evaluation. Exceptions are JW 681, JW 876, and JW 959 
where no photometric extinctions could be measured and our 
best estimates for Ay come from spectroscopy with an uncer- 
tainty determined from the fitting error ellipse that is similar to 
the uncertainty in k. 

3.4. Luminosity, effective temperature, and radius 

Luminosities L„ of our target components were derived from 
bolometric magnitu des by applying bolometric corrections BCj 
dHartigan et al.lll994l) to the measured /-band magnitudes with a 



10 



Daemgen, Correia, & Petr-Gotzens: Protoplanetary disks of T Tauri binaries in the ONC. 



-2 





1 1 1 1 


i i 1 1 


1 1 1 i i i I i i 1 1 1 i i i 

2-0OM Q 


i i 


■ 






-' \l .40M„ 


■ 






■ 






: 


- 




■ 






- 


■ 










- 


': 








■to 




























'■_ 






















A\ \o.10M 


i 










." \\ - O.IMyr 


: 

: 








v\\n 


'0 3 My- " 

MfjfT 




I 


variability 






II 0.0 My r 












v 30.0Myr 






h random error 










. i , , 


i i . . 




. i 


lOO.OMyr 
1 1 ■ ■" 



3.9 3.8 3.7 3.6 3.5 3.4 

iog(T eff ) 

Fig. 6. HR-diagram with evolutionary tracks from Sies s et al.l 
(2000). Open symbols are target components with high veil- 
ing r% > 0.2, filled symbols have r# < 0.2. Typical uncer- 
tainties from random errors and the resu lting uncertainty from 
0.2 mag intrinsic photometric variability (Carpe nter et al.ll2~00ll) 
are shown in the lower left. For the determination of the parame- 
ters (age & mass, see Table |6), more tracks and isochrones were 
used than shown in the plot; they were omitted for a clearer il- 
lustration. 

distance to Orion of 414 ± 7 pc dMenten et al.l|2007l) and /-band 
extinctions of Aj — 0.27Ay (ICohen et al.l 198 lh . The /-band was 
chosen to help us minimize the impact of hot circumstellar mate- 
rial, which mainly contributes flux at longer wavelengths (K and 
L-band) i.e. closer to the maximum of the T ~ 1500Kb lack- 
body emission from the inner dust rim (Me yer et al.ll 19971) . The 
resulting luminosities are listed in Table [6] along with effective 
temperatures r e (j derived from their spectral typ es and SpT-/ e ff 
relati ons (earlier than M0 : [Schmidt- Kalerll982t later or equal to 
M0: lLuhman et al.1l2003l) . Luminosity uncertainties were propa- 
gated from the magnitude, extinction, and distance uncertainties 
but do not include any of the intrinsic variability of the targets 
since we cannot estimate the magnitude of the effect for any in- 
dividual target componenfl We estimate, howe ver, an average 
impac t of variability by varying Mj by 0.2 mag ( Car penter et al.l 
1200 lh and rederiving the luminosity of each target. We observe a 
difference in the derived luminosities of up to 20% with a sample 
median of 0.08 dex. This is twice as large as the typical propa- 
gated random uncertainties. Stellar radii R were then calculated 
from L* and / e ff. 

3.4.1 . The HR diagram: Ages and masses 

Effective temperatures and luminosities were used to derive 
ages and masses by comparison with evolutionary tracks from 
Siess et alj(l2000l) . Fig.[6]shows the position on the HR-diagram 
of all components with measured T s g and L„. The derived 
masses are - except for three targets - below 1 M Q , and ages 
are found to be in the range of 10 4 -10 7 yr, with an average of 
0.97 Myr. One object, the secondary component of [HC2000] 73, 

2 The variabi l ity of a few of the target binaries was observed by 
Carpe nter et all ( 120011) . although they did not resolve the individual 
components. 



has an estimated mass of 0.09 ± 0.05 M Q , indicating that it is a 
possible substellar object. Th e details of this binary are discus sed 
in an accompanying paper dPetr-Gotzens & DaemgenlfeOl ll in 
prep.). 

It is apparent that the targets that are classified as the 
youngest, are those with the highest veiling. Their high lumi- 
nosities thus probably do not represent extreme youth but are 
rather an indication of hot circumstellar material contributing 
near-infrared flux even in /-band. This is not properly accounted 
for when extracting ages from the HR diagram because the 
infrared-excess in the /-band is unknown for our targets and 
therefore not subtracted from their brightness. A consequence is 
the apparent non-coevality of binary stars with at least one com- 
ponent with high veiling, as discussed in Sect. 14.1.21 and Fig.fTOl 
Furthermore, apparently old ages can be caused by underesti- 
mated extinctions, making these targets appear underluminous 
and thus too old. However, since the evolutionary tracks for stars 
of a certain mass are almost vertical (T e ff » const) in this part of 
the HR diagram, the uncertainty in luminosity does not translate 
into an equally large uncertainty in mass. 

3.5. Accretion: W(Bry), L acc , and M acc 

The accretion activity of each target component was inferred 
from the Bry emission feature at 21665 A. Our measurements 
are described in this section. 

3.5.1 . Equivalent widths of Bry emission 

We measured the equivalent width of the Bry line 

W Bry = f ^^dA (4) 

JBvy r c 

in an interval of width 56 A around the center of Bry (which was 
individually fit and recentered between 2 1654 A and 2 1674 A) 
for all target components. The integration interval was chosen to 
ensure that we include as much line flux as possible, while min- 
imizing the influence from the continuum noise around the line 
at the spectral resolution of our observations. The continuum F c 
was determined from a linear fit to the local pseudo-continuum 
in a region of width 130 A shortward and longward of the inte- 
gration limits. The above definition of equivalent width (eq. [4]) 
returns positive values for emission lines and is negative in case 
of absorption. 

To assess the significance of Bry emission, we inferred the 
noise level at the position of Bry from the surrounding contin- 
uum. This allowed us to derive a measure for the probability of 
the presence of a gaseous accretion disk in a binary component 
(see AppendixlAli. 

The procedure of measuring equivalent widths was per- 
formed in two stages. In the first stage, the method was applied 
to the uncorrected spectra to serve as a measure of accretion 
disks (Appendix lAli. where we sought to assess the significance 
of the measured emission features and thus do not wish to in- 
troduce additional uncertainties by means of further modifica- 
tions of the spectra. However, spectral features were reduced in 
strength when veiling was imposed on the spectrum and thus 
the Bry equivalent width had to be corrected for r% to represent 
the actual emission emitted from the accretion process. This re- 
quired a second step of either calculating Wa = W^ asured -(1 +r^), 
as inferred from eqs. (Q]i and d4), or applying eq. {TJ with r# val- 
ues from Table [6] to the reduced spectra (for which extinction 
does not change the equivalent width and can be assumed to be 
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Fig. 7. Veiling-corrected Bry equivalent width as a function of 
NIR excess in (H-K s ) color. Circles represent primary, and 
diamonds secondary binary components. The detection limits 
of target components with insignificant equivalent widths are 
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equal to zero for this calculation). Both options returned very 
similar results. We chose to modify the spectra and remeasure 
V^Bry to also obtain good estimates for the continuum flux noise, 
which does not necessarily transform in the same way. The re- 
sults for this measurement of the actual equivalent widths of the 
Bry emission are listed in Table [6] When r# was unknown, no 
correction was applied and values in Table[6]correspond to lower 
limits of Wsry. 

3.5.2. Bry equivalent widths versus NIR excess 

We compared the color excesses Eh-k s , which measure the ex- 
isten ce of hot circum stellar material around each binary compo- 
nent (ICieza et al.l2005l) . with our veiling-corrected Bry emission 
values. Fig.|7]shows that, while all but one target with significant 
emission in Bry exhibit a NIR excess, the opposite is not true and 
many targets with NIR excess are found that do not show signif- 
icant signs of hydrogen emission. This imbalance is well-known 
and discussed further in Sect. I4.2. II 

The figure shows some interesting features. It seems that 
there is a slight systematic offset in the calculated excesses rela- 
tive to the origin, as we see a number of targets clustering around 
Eh-k s ~ -0.03. Assuming that these are targets with no sig- 
nificant color excess, we can infer a small mismatch between 
the color scales of the theoretical and measured values used to 
derive Eh-k^- This is partly a consequence of using dwarfs for 
the theoretical colors {H-K s ) t h eor instead of pre-main sequence 
stars. W hen using the colors of pre-main sequence stars de- 
rived by iLuhman et alJ (l2010h . we measured an average shift 
of 0.015 mag towards redder H-K s colors that moved the ac- 
cumulation of low color-excess compon ents towards th e origin. 
Since this correction was small and the ILuhman et al.l spectral 
sequence does not cover all spectral types of our sample, we use 
color-excesses derived from the dwarf colors. In addition, the an- 
ticipated systematic offset was negligible compar ed to the possi- 
ble ph otometric variability of ~0.2 magnitudes ( Carp enter et al.l 
1200 ll) . Although variability does not change our (qualititative) 



conclusions, it might be a reason for the position of the targets 
with negative excess with and without significant Bry emission. 

Targets in the bottom right quadrant display Bry in absorp- 
tion; since we also detect a NIR emission excess, there is a possi- 
bility that Bry emission generated by accretion cancels with part 
of the absorption feature. To estimate the real emission strength, 
one could measure and subtract the strength of Bry absorption 
from photospheric standards of the same spectral type. This con- 
cerns, however, only the components of JW 260, a binary that 
was excluded from the discussion and conclusions owing to its 
earlier spectral type. We thus skip a more thorough evaluation of 
the accretion state of this binary here. 

3.5.3. Bry line luminosity and mass accretion rates 

For target components with full knowledge of extinction and 
veiling, we calculated Bry line luminosities and the mass ac- 
cretion rate M acc . The accretion luminosity was derived through 

log(L acc ) = (1.26 ± O.19)log(L Bry /L ) + (4.43 ± 0.79) (5) 

dMuzerolle et aljH998l) where the Bry line luminosity is defined 



: Anr 



JBry 



(F A -F c )dA 



(6) 



in the same integration limits used for the equival ent width and 
using a distance to the ONC of r — 414 + 7pc dMenten et al.l 
2007). To measure the luminosity of the Bry line, the spectra 
must be flux-calibrated. This was achieved by comparing the 
A^-band photometry from Table [5] with synthetic photometry 
obtained by convolving our measured spectra with the 2MASS 
K s filter curve and integrating with zp{F A ) = 4.283 x 10~ 7 for 
the zeropoinfl Since the filter curve extends to slightly bluer 
wavelengths than our NACO spectra, we had to extrapolate the 
spectrum. To estimate the impact of our linear extrapolation on 
the integration result, we assumed several extrapolations with 
different slopes (up to unreasonable values). The resulting vari- 
ation in LBry is small, mainly because the extrapolated part co- 
incides with the steep edge of the filter and was kept as an ad- 
ditional uncertainty. After correcting the calibrated spectra for 
veiling and extinction, the data was multiplied with the filter 
curve, integrated over Zf s -band, and converted to L acc according 
to eq. (0. Uncertainties were estimated from error propagation 
of all involved parameters including the extrapolation error and 
the empirical uncertainties in eq. A source of additional un- 
certainty that we could not quantify in greater detail from our 
observations is variability. Since the photometry was not taken 
simultaneously with our spectral observations, the calibration of 
our spectra might suffer from additional uncertainty when these 
young targets were observed in different states of activity. Since 
we had no estimate of the size of this effect for an individual tar- 
get, we did not introduce any correction but wish to caution that 
individual accretion luminosities might be offset from the true 
value. However, assuming that the effect of variability in the ac- 
cretion luminosities is random, the sample statistics should not 
be biased. 

The mass accret ion rate was calculated according to 
iGullbring et all (1 19981) as 



^acc ^* 



GM, \R in -R, 



(7) 



http://www.ipac.caltech.edu/2mass/releases/allsky/doc/sec6_4a.html 
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for a stellar radius R* and mass M, from Table [6] and the gravia- 
tional constant G, assuming that material falls onto the star from 
the inner rim of the disk at R[ n ss 5R*. The accretion luminosities 
and mass accretion rates are listed in Table [6] 



3.5.4. Is the Bry emission generated by magnotespheric 
accretion? 

The main source of Bry emission in T Taur i stars is often as - 
sumed to be magnetospheric accretion (e.g. iBeck et alJ l2010). 
However, mechanisms such as stellar wind, disk wind, outflow, 
or photoevaporation of the disk by a nearby high-mass star 
can al so contribute to the Bry emission observed in low-mass 
stars (lHartmann et al.1 Il990t lHartigan et"ai] 1 1995b lEisner et all 
120101) . While the emission region from magnetospheric accre- 
tion should be located close to the stellar surface, the sources 
of most other mechanisms are expected to be located further 
away at several stellar radii or even in the outer parts of the disk. 
Detecting a spatial displacement of the Bry emitting region from 
the star locus would refute the possibility of magnetospheric ac- 
cretion as the origin of the emission and render other explana- 
tions more likely. 

To test this, we used the sky-subtracted raw frames of the 
spectral observations, showing the spectral traces of both com- 
ponents of each target binary. The orientation of the images is 
such that the dispersion direction is roughly aligned along the 
columns of the detector, while the spatial information is aligned 
along the rows. In each row, we fit two Gaussians to the two pro- 
files of the component spectra, thus measuring the spatial loca- 
tion and width of the spectral profile of the targets in each wave- 
length bin of size ~5 A. Depending on the S/N of the individual 
observations, the locations of the trace center can be determined 
to an accuracy of 0.008-0.10 pixels, which, at the spatial reso- 
lution of the observations of ~0.027 arcsec/pix, corresponds to 
~0. 1-1 . 1 AU at the distance of the ONC. 

We found no significant offset at /l(Bry) from the rest of the 
trace for any of the targets in which we observed Bry in emis- 
sion. Neither did we detect extended emission in excess of the 
width of the spectral trace at similar wavelengths. This indicates 
that the emission indeed originates from a small region close 
to the stellar surface and is likely to be the product of magneto- 
spheric accretion. However, extended or displaced emission may 
well be generated mainly perpendicular to the slit. In this case, 
a displacement of the emission peak could be observed in the 
dispersion direction. Since we did not measure frequency shifts 
of spectral features, we cannot exclude this possibility for any 
individual target component. 

While the Bry emission seems to come from a region close 
to the star for all targets, we observed one target (TCC55) 
in which emission at /KH2) = 21218 A comes from an ex- 
tended region around the star, rather than the star itself. We ob- 
served at least three diffuse H2 emitting regions along the slit, 
each several pixels wide, one of them apparently surrounding 
the binary. We used K-b and images from the MAD instrument 
dPetr-Gotzens et alJ2008l) to investigate the surrounding area for 
possible emission sources close to this target located only 0.'49 
from 9 1 OriC. We identified bow-shocks coinciding with the lo- 
cation of the H2 features in the spectra. The diffuse emission 
around the binary itself might represent the remaining material 
of a proto-stellar envelope, which would be indicative of a young 
evolutionary state (i.e. Class I) of the binary. 
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Fig. 8. Number of target components per spectral type. Gray ar- 
eas: primaries. Black outline: secondaries. The distributions of 
primaries and secondaries are only slightly different in spectral 
types. 

4. Discussion 

4. 1 . Stellar parameters and sample biases 

To detect possible biases in the numbers derived from the com- 
ponents of our binary sample, we discuss the degree to which 
the populations of primaries and secondaries differ and whether 
they are typical of the ONC population. 

4.1.1. Spectral types 

It is known that the strengths of accretion parameters correlate 
with the mass of a star, which (considering the limited spread in 
age) can be represented by the typically more tightly constrained 
spectral typqj In Fig. [8] we see that the spectral type distribu- 
tion of the primaries peaks at slightly earlier spectral types of 
M0-M1 than the secondaries (M2-M3). This difference is sig- 
nificant at 98%, according to a Kolmogorov-Smirnov (K-S) test. 
However, both distributions are indistin guishable from the spec- 
tral types of the entire ONC population ( Hill enbrand! 1 997 ). with 
K-S probabilities for different distributions of 65% and 57% for 
primaries and secondaries, respectively. This means that both 
primaries and secondaries are 'typical' members of the ONC, 
whereas the primary and secondary spectral type distributions 
deviate slightly and differences in the derived parameters (such 
as accretion rates and Bry-emission strength) can partly be at- 
tributed to the - on average - earlier spectral types of the pri- 
maries. 

4.1.2. Relative extinctions and ages 

Interstellar extinction through embedding in the Orion molecular 
cloud expresses itself as a spatially variable source of extinction 
that is, nevertheless, very similar for all components of a stellar 
multiple. An additional source of extinction that can be very dif- 
ferent even for components of the same binary can be caused by 

4 That accretion sig natures do correlate with spectral type can be seen 
from e.g. Table 7 in lWhite & Ghezl d200lh . 
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Fig. 9. Extinction of the primary versus versus secondary com- 
ponents. Circles show targets with Ay determined from the 
dereddening of the CTT locus. Squares show the spectroscop- 
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because photometric extinctions could not be measured. Filled 
symbols indicate targets with both components having a low 
veiling rx < 0.2. The dashed line corresponds to equal extinc- 
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Fig. 10. Relative ages for the primary and secondary components 
of all target binaries where both components could be placed in 
the HR diagram. Filled squares show targets where both com- 
ponents exhibit low veiling r# < 0.2. Targets with open squares 
have at least one component with strong veiling > 0.2, thus it 
is likely that the age estimation from the HR diagram is biased, 
since extra luminosity from accretion makes the targets appear 
brighter and thus younger. 



circumstellar material such as a (nearly edge-on) circumstellar 
disk, obscuration by the other component's disk, or a remaining 
dust envelope. 

In Fig. [9] we can identify and verify the impact of the dif- 
ferent sources of extinction. Binaries composed of two low- 
veiling components have similar extinctions, according to the 
level of embedding in the cloud, whereas binaries with high- 
veiling components do not display such correlation. This might 
be due to dust material from a disk that is detected by means of 
its hot-conti nuum emissio n from magnetospheric accretion, i.e. 
veiling (e.g. lBouvier et alJl2007h . In particular, the magnitude of 
extinction is a strong function of the angle under which the sys- 
tem is observed with close to edge-on disks causi ng a strong red- 
dening of the NIR colors dRobitaille et alJ l2006). However, our 
observations do not enable us to determine to the inclinations of 
the circumstellar disks, hence conclusions about disk orientation 
and alignment cannot be drawn. 

The stellar components within a binary are close to being 
coeval. Fig. [10] shows a well-defined correlation of primary and 
secondary ages for those binaries with little or no veiling. The 
five binaries with at least one component of > 0.2 are clearly 
located off the sequence of coeval binaries, which is probably 
due to a non-negligible amount of veiling in the /-band for tar- 
gets with high rK values. In our present study, we did not attempt 
to derive accurate absolute ages for our sample, but testing for 
equal ages within binaries does serve two purposes: i) A san- 
ity check. Assumin g that binary components d o form reason- 
ably close in time (Kraus & Hill enbrand! 12009). we confirmed 
that our derived r e (j and L* were accurately determined because 
they result in reasonable values when placed in an HR-diagram. 
We inferred that the derived parameters (e.g. L acc , M acc ) are of 
sufficient quality to help us derive the conclusions of this paper. 
ii) The derived parameters have no age dependence. Since there 
is no systematic difference between the primary and secondary 
age of the binaries, we were able to exclude any dependence on 



age of the derived relative parameters between the populations 
of the primaries and secondaries. Hi) Coevality is consistent with 
the physical binarity. 

While we derived ages r = log(age) that span a range of 
~ 5.5 < t <~ 6.5 even for the well-behaved class of low-r# bi- 
naries, this can pro bably not be attributed to a real age spread. 
I Jeffries et alJ (|201 1 ) observed no difference in age between stars 
with and without disks in the ONC, indicating that the age spread 
must be shorter than the disk lifetime, i.e. significantly more con- 
fined than traditionally assumed. We conclude that not only ob- 
servational uncertainties, but also an intrinsic scatter in the lu- 
minosities at constant age must be present. Despite the impor- 
tance of this observation for the absolute ages of members of 
star forming regions, we assumed that our relative age measure- 
ments apply, since the boundary conditions are very similar for 
both components of the same binary. 

4.2. Disk evolution around the components of visual ONC 
binaries 

The absolute and differential abundances of disk signatures in 
the components of our target binaries were found to be a function 
of binary parameters, as we discuss in the following sections. 

4.2.1 . The disk fraction of binary components 

We derived the fraction of ONC binary components harboring 
an accretion disk. This number was compared to the single-star 
disk frequency in the ONC and to binary samples of other star 
forming regions to expose the effect of binarity and cluster envi- 
ronment on the evolution of circumstellar disks. 

The probability density function of the disk frequency was 
derived in an Bayesian approach as well as the probability of 
an individual binary component harboring an accretion disk 
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(see Appendix lAl individual disk probabilities in Tabled. As 
the strength of disk signatures depends on spectral type (e.g. 
IWhite & Ghezll2001l) . we excluded the binary JW 260 from fur- 
ther evaluation of disk frequencies, since its spectral type is con- 
siderably earlier than for the rest of the sample. We measured 
an accretion disk fraction of F = 35+g% in our sample of 42 
spectroscopically observed binary components ( Appendix lAb. 

Analyses of the Bry line are known to return a smaller frac- 
tion of accretors than the Ha emission feature, which is typi- 
cally used to decide whethe r an individual star can be classi- 
fied as a classical T Tauri star. Folh a & Emersonl fcOOl) observed 
a number of classical T Tauri stars in Taurus-Auriga, measur- 
ing NIR emission-line strengths including Brackett-y signatures. 
Twenty-four of their targets are in the range of K3-M6 spectral 
types, comparable to our sample, and three show no emission in 
Bry. Hence, a fraction of / = 0.125 of all classical T Tauri tar- 
gets would have to have been misclassified as weak-line T Tauri 
stars from their Bry emission. We thus expected a number of 
//(I _ /) x (F ■ AO ~ 2 (with N - 42 target components) classi- 
cal T Tauri stars to be classified as non-accreting. The corrected 
fraction of classical TTauris among binary star components in 
the ONC was thus Fctt = 40^„°%. Since many studies refer 
to the fraction of Ha-detected classical T Tauri stars rather than 
the number of accretors from Bry emission, we used the latter 
numb er to compare it with st udies of accretion disk fre quencies. 

IHillenbrand et alJ l l 19981) and iFuresz et all (120081) found a 
frequency of accretion disks bearing single stars in the ONC 
of 50% and 55%, respectively. Both are marginally (lcr and 
1.5<x) larger than our measured fraction of 40^°%. While the 
IHillenbrand et al.l sample was de rived from th e Ic — K color in- 
stead of Ha measurements, the Fur esz et al.l sample is biased 
towards classical T Tauri stars making their estimate an up- 
per limit of disk frequency. Both of these findings prevent us 
from drawing firm conclusions about the difference between the 
disk frequencies of single stars and binaries. This evidence of 
a lower disk frequency around 100-400 AU binary components 
will hence need future confirmation from observations using 
comparable diagnostics (preferably Bry) in an unbiased com- 
parison sample. 

Similarly, we found evidence of an underrepresentation of 
dust disks in Orion binaries. Sixteen out of 27 target compo- 
nents (excluding JW260) with measured H-K s colors show 
signs of a dust excess, that is a fraction of 59+15% compared 
to the values of 55%-90 % found for single stars in the ONC b y 
IHillenbrand eTail (Il998l) and 80+8% from lLada & Ladal (120031) . 
Since these numbers were derived using indicators other than the 
H-K s excess, which is known to typically return a comparably 
small fraction of dust disks compared to e.g. K-L (IHillenbrand! 
120051) . this evidence cannot be quantified in greater detail. 

Our numbers suggest that there is a higher frequency of tar- 
get components with inner dust disks (59+15%) than accretion 
disks (40*1°%). This discrepancy agrees with ob servations of 
single stars in various star-forming regions, where iFedele et al.1 
(2010) concluded that accretion disks decay more rapidly than 
dust disks in a particular cluster. 

The presence of dust and accretion disks around the binary 
components of our sample is thus consistent with i) observations 
of single stars in a variety of young clusters that dust disks are 
more abundant than accretion disks in the same cluster and ii) 
the expectation that disk lifetimes are shorter for disks in binary 
systems than for single stars with comparable properties. The 
latter is also theoretically mo tivated by the missing outer disk 
through dynamical truncation dArtvmowicz & Lubowl 1994b and 
the resulting reduced feeding of the inner disk from the outer 
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Fig. 11. Histograms of binary separation as a function of com- 
ponent accretion-type (the y-axis tickmarks indicate one binary 
each). Most of the close binaries are of type WW or CC, i.e. 
synchronized in their disk evolutionary state. 

disk material (Monin et aT1 l2007l) . As expected, disk frequency 
is a fu nction of binary se paration (also see Sect. I4.2.3I I. For ex- 
ample, ICieza et ail d2009l) found the component disk frequency 
in tight <40AU binaries of Taurus-Auriga to be significantly 
lower, at less than one half of the single-star disk frequency. Both 
their wider binaries (40-400 AU) and our sample (100-400 AU), 
however, have disk frequencies lower but comparable to single 
stars, indicating that there is a separation dependent mechanism. 

4.2.2. Synchronized disk evolution in ONC binaries 

We detected a significant overabundance of close, <200 AU pairs 
of equal emission state, i.e. with both components accreting 
(CCjfl or both components showing no accretion (WW) signa- 
tures. This is apparent in Fig. inland a K-S test indicates 99.5% 
probability that the separation distribution of equal pairs (CC 
and WW) differs from the separations of mixed pairs (CW and 
WC). 

To investigate a possible correlation between the synchrony 
of disks in ONC binaries and their separations, we split our sam- 
ple into binaries with projected separations larger and smaller 
than 200 AU. From the content of accreting and non-accreting 
components in the two separation bins, we predicted the average 
number of CC, WW, and mixed systems by random pairing, and 
compare it to our measured distribution. If there were no corre- 
lation between the evolution of both components of the same bi- 
nary, the randomly paired sample of the same number of compo- 
nents (W and C) should be consistent with our observed sample. 
With 18 non-accreting components and 6 accretors in the sample 



5 To ease the reading of the binary categories, we use the common ab- 
breviations 'C for accreting components (referring to classical T Tauri 
Stars), and 'W for non-accreting components (weak-line T Tauri stars) 
where in the designation of a binary pair (e.g. 'CW') the first position 
describes the primary (here 'C') and the second the secondary (here 
' W) component state. We note, however, that we only refer to the pres- 
ence of accretion as measured through Bry, which is correlated with but 
not equal to the the distinction between weak-line and classical T Tauri 
stars (see also the discussion Sect. l4.2.Tt . 
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of 12 binaries with separations <200 AU, we expect an average 
of ~1.2xCC, ~6.9xWW, and ~4.1 mixed systems, as predicted 
by a Monte Carlo simulation. However, we found 3+lxCC, 
9+0xWW, and 0*2 mixed pairfl This is clearly incompatible 
with the prediction of random pairing. On the contrary, wide- 
separation binaries are well described through random pairing 
with predicted values of 2.2xCC, 4.1xWW and 4.8xmixed and 



measured values of l^xC C, 3^xWW and T^xmixed. 

IWhite & Ghezl d200lh observed a similar underdensity of 
mixed pairs among binaries with separations <210 AU in a sam- 
ple of 46 binaries in the Taurus-Auriga star-forming associa- 
tion. They concluded that synchronized evolution, which they at- 
tributed to the existence of a circumbinary reservoir, can more or 
less equally replenish the circum primary and circumse condary 
disks, as previously suggested by Prato & Simon] d!997l) . 

Can circumbinary disks of sufficient size survive in the ONC 
and thus be the cause for the synchronization of disk evolu- 
tion? The typi cal size of a disk in th e Trapezium region is be- 
low 200 AU (IVicente & Alvesl 120051) and only 3 of 149 ana- 
lyzed systems - the authors claim completeness for large disks 
> 1 50 AU at moderate extinctions - were found to have disk sizes 
>400AU. However, circumbinar y disks have inner radii of a t 
least twice the binary separation (Artvmowicz & Lubowl [l994). 
This requires circumbinary disk sizes of more than 400 AU in 
diameter for a binary with 100 AU separation and even >800 AU 
circumbinary disks for 200 AU binaries. Assuming that dy- 
namical interaction s ( Olczak et alj|2006h and photoevaporation 
(M ann & Williams! 120091) are the reason for disk truncation in 
the ONC, the observed size limits of single star disks should also 
apply to circumbinary matter. This implies that circumbinary 
disks of binary systems >100 AU should not be largely abun- 
dant since they are typically truncated to radii below the dynam- 
ically induced inner hole radius. These considerations render it 
unlikely that stable circumbinary disks are the reason for disk- 
synchronization in the inner regions of the Orion Nebula cluster. 

In agreement with this, none of these systems th at were ob- 
serve d with circ umbinary material (e .g. GG Tau, Dutr evet alj 
1994 V892Tau. iMonnier et all 120081: Orion proplyd 124-132, 
Robberto et al.l 120081) have binary separations of more than 



100AU. This points to the possible universality of the result. 
IVicente & Alvesl J2005 ) observed no trend in the disk sizes with 
distance to 6 1 OriC out to ~4', which means that Orion disks 
are small, independent of their position in the inner ~50 arcmin 2 
of the cluster. Furthermore, we would expect to see a larger 
ratio of CC binaries to WW binaries at large distances from 
9 l Ori C when circumbinary disks are present only in the outer 
parts of the cluster and transfer a significant amount of material 
to the individual stellar disks. The distribution of CC and WW in 
our ONC does not, however, increase with distance to 6 l Ori C 
(Fig. FPU) . Finally, although in Taurus dAndrews & Williams! 
120071) more disks with large radii have been observed than in the 
ONC, probably owing to its weaker dynamical interactions and 
irradiation, the similarity of the observed parameters (such as the 
200 AU limit for synchronization) seems to suggest that the disk 
feeding mechanism in Taurus binaries is probably similar to that 
in Orion, i.e. not due to replenishment from circumbinary disks. 

What other mechanism could synchronize the circum- 
primary and circumsecondary disks in <200AU systems? 
Since mass accretion rates are a function of stellar mass 
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Fig. 12. Histograms of distance to 9 l Ori C as a function of ac- 
cretion type. The histogram shows no indication that binaries of 
any type are more likely at any distance. The v ertical dashed 
line shows the radius (460"; Rei purfh et al.ll200"7b inside which 
the ratio of wide binaries (075- T/5) to close binaries (07 15-075) 
drops considerably, probably owing to dynamical interaction. 



6 uncertainties are derived from the possibility of a binary changing 
classificiation within its Icr limit of Wb^,, i.e. a CC might turn into a 
CW, if its secondary is classified as C but with an equivalent width that 
is less than lcr away from W ra ; n . 



(IWhite& Ghezl 1200 1± and disk truncation radii are similar in 
equal mass systems (lArtvmowicz & Lubowl 19941) . synchroniza- 
tion of disk evolution might arise, if close binary systems are 
preferentially equal mass systems. For 13 binaries of the sam- 
ple, we were able to derive masses of both binary components 
(see Fig. [T3l . For these, we observe that all (four) systems with 
mass ratios of q>0.8 and separations <200 AU are of WW type, 
agreeing with the hypothesis of high mass ratios being the cause 
of synchronized disk evolution. The significance of this result, 
however, is low. There is a 21% chance that the mass ratios of 
close (<200AU) systems are drawn from the same parent dis- 
tribution as mass ratios of wider pairs (K-S test). A larger sam- 
ple of spatially resolved spectroscopic observations of pre-main 
sequence binaries is needed to decide whether mass ratios are 
the main driver of the synchronization of binaries closer than 
200 AU. 



4.2.3. Differential disk evolution in binaries 

Mixed pairs with accreting (CW) and non-accreting primaries 
(WC) - within the uncertainties - are equally abundant: CW 
pairs appear 4+2 times while WC are measured 3+1 times. 
Although not statistically significant, this is evidence against a 
strong preference for primaries to have longer lived disks than 
the less massive secondary. However, longer lived disks around 
primaries are suggested by the ory since disks around second aries 
are truncated to smaller radii (lArtvmowicz & Lubowll 19941) and 
dissipation times are predicted to scale like R 2 ~" with R the disk 
radius and q s al-1.5 dMonin et al.l l2007, and references therein). 
iMonin et alJ (120071) found that their measured overabundance of 
14xCW versus 6xWC is consistent with this effect taking place, 
however, with other factors (i.e. initial disk conditions) having 
a more a dominant impact on the lifetimes of circumprimary 
and circumsecondary disks than the differential scaling with R, 
which is only strong for binaries with low mass ratios q<0.5. 
Our data agree with this proposed weak correlation of the bi- 
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Fig. 13. Mass ratios of binaries as a function of projected binary 
separation, indicating whether both binary components are ac- 
creting (filled circles, CC), neither component is accreting (open 
circle, WW), or either component shows signs of accretion (half- 
open circles = CW and WC). This plot only contains 13 targets, 
since masses could not be derived for all binary components in 
the sample. For some target binaries with mass ratios close to 
1, the primary (as estimated from the NIR colors and listed in 
Table |5]l turned out to be the less massive component. The mass 
ratios of those targets were calculated as the inverse leading to 
^-values <1. 



nary mass ratios with the abundance of CW-binaries in Orion, 
though this result is limited in significance by the small number 
of mixed systems in our sample. 

It is noteworthy that the existence of mixed pairs, together 
with their property of having wider separations (Fig.fTTTi. can in- 
troduce difficulties in the interpretation of binary stu dies that do 
not re solve their targets into separate components. ICieza et all 
(2009) used NIR photometry of unresolved binaries with known 
separations from several star-forming regions (Taurus, p-Oph, 
Cha I, and Corona Australis) to infer a smaller separation in bi- 
naries with no accreting components than in accreting binaries. 
Besides the proposed shorter disk lifetimes around close binary 
components, there is an alternative interpretation of their data 
that they did not discuss. Since they did not resolve binaries into 
separate components, they were unable to distinguish CC, CW, 
and WC-type binaries, but merged them all into the category of 
having at least one disk. As we showed earlier, however, the sep- 
aration distribution of CW and WC binaries differ significantly 
from equal-accretion binaries including CC (Fig. fTTT i. When 
joining the three categories with at least one accreting compo- 
nent, the combined separations display a distribution with on av- 
erage larger separations than the WW distribution. Only from 
our resolved population, can we interpret this as a lack of mixed 
(CW+WC) pairs and not a lack of accreting (CW+WC+CC) 
components in close binaries. 

4.2.4. Accretion luminosities and mass accretion rates 

Figs. Q4] and Q3] show a luminosity histogram and component 
mass-accretion rates as a function of stellar mass, respectively. 
In binaries with two accreting components, it is usually the 
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Fig. 14. Histogram of the accretion luminosities of the primaries 
(gray shaded area) and secondaries (hatched) calculated from the 
Bry emission. For comparison, the accretion lumin osity distribu- 
tions of sing le stars of Orion fro m lRobberto et al.l d2004 dashed 
outline) and D a Rio et al.l d2010l dotted outline, scaled by a fac- 
tor of 1/20 for clear comparability) are overplotted, both lim- 
ited to the same range of stellar masses as in our binary survey. 
The distributions are slightly offset relative to each other to make 
them more visible. 




Fig. 15. Mass accretion versus stellar mass for all significant 
emitters of the sample (filled symbols) and upper limits to 
all other targets with measured M acc and M* (open symbols). 
Primaries are marked with circles, secondaries with diamonds. 
Asterisks show the m ass accretion rates of single stars in Orion 
(Robberto et alj|2004h . where as plus signs and up per limits are 
binary components in Taurus (I White & Ghedl2001l) . 



more massive component that has the higher accretion luminos- 
ity. Accordingly, we see a tendency for the subsample of pri- 
mary components to have slightly higher relative accretion lu- 
minosities L acc /L* than the secondaries of our Orion binaries. 
The derived accret ion luminosities agree very well with s ingle 
stars in the ONC dRobberto et all 12004 iDa Rio et al]l2010h for 
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log(L acc /L ) > -1.5, which is the average sensitivity limit of the 
accretion luminosities derived from our Bry observations. 

Similarly, the mass accretion rates of the ONC T Tauri bi- 
nary components as a function of stellar mass are (except for 
three outliers) comparable to those of single stars in the ONC. 
Fig. [15] shows that Orion single stars occupy almost the same 
area in the log(M acc )-log(M») diagram, with tendency towards 
slightly lower accretion rates. This tendency i s most likely an 
observational bias, since Robbertp_etalJ(2004) use [/-band ob- 
servations with the Hubble Space Telescope that are more sen- 
sitive to lower mass accretion rates than our Br y data. We also 
overplot stellar components of Taurus binaries ( Wh ite & Ghezl 
1200 ll) . which have similar mass accretion rates. 

Comparable mass accretion rates of singles and binaries in 
the ONC and Taurufl are not self-evident, considering the differ- 
ent disk populations of the ONC and Taurus. Since disk masses 
in binaries are lower than in single stars of the same star-forming 
region because of disk truncation ( Artvmo wicz & Lubowf l994) 
and disks in Orion are less massive than disks in Taurus 
(M ann & Williams 2009), uniform accretion rates indicate either 
different disk lifetimes or variable efficiency for the replenish- 
ment of an existing disk. The hypothesis of shorter disk lifetimes 
would corroborate the evidence from Sect. l4.2.2l that we observe 
fewer disks around binary components than were measured for 
single stars. 

Three binary components (TCC52A&B, JW391A) have 
comparably high mass accretion rates to the rest of the sample 
in Fig. [15] The responsible mechanism is obscure, since the two 
respective binaries have no remarkable properties in common 
relative to other targets. Their mass, separation, and distance to 
9 l OriC are unremarkable. While JW391A shows no peculiar- 
ity in luminosity, the two components of TCC 52 are found to be 
the most luminous targets with respect to their mass. The derived 
very young ages and large radii might indicate an earlier evolu- 
tionary stage, i.e. lass I, which would agree with the higher mass 
accretion rates of TCC 52A+B than to older class II components 
(Robitaille et al.l2006b . Regardless of their properties, it remains 
possible that all three components were observed in a temporary 
state of high activity. 

To assess whether disks in binaries need significant replen- 
ishment to survive in sufficient quantities, we estimated disk life- 
times = Mdi s k/M acc from the ratio of disk mass to mass 
accretion rate and compare it to the age of the star-forming re- 
gion. To estimate M^sk, we compiled upper limits to the total 
mass of dusty material around both binary components from the 
literatur e of millimeter observatio ns (JW519, JW68 1, TCC 52 
TCC 97. lMann & WilliamsbOlOl: TCC 52, TCC55. lEisner et all 
l2008h . 

Only for TCC 52 do we have available both mass accre- 
tion rates and the total mass of the surrounding (disk) ma- 
terial. The two measurements of total disk mass of TCC 52 
disagree at the l c r-leve l: 0.0288+0.0029 M Q are derived by 
Man n & Williamsl (1201 Oh and 0.042+0.009 M by lEisner et all 
(2008). Nevertheless, we now illustrate that an "order-of- 
magnitude" estimate is possible when we assume a total 
disk mass of ~0.03M o . Disk radii can be estimated from 
their dynamical truncation radii to be ~0.38 and ~0.3 times 
the binary separation for primary and secondary, respectively 
(lArmitage et al.ll 19991) . considering ou r derived binary mass ra - 
tio of q * 0.65 and M disk oc R disk (e.g. lMann & WilliamslbOlOh . 
This results in individual disk masses of M^j™ « 0.017 M Q and 



.Mdisk ~ 0-013 M Q . The derived disk lifetimes for these two target 
components with the highest mass accretion rates of our sample 
are t^™ * 1 .4 X 1 4 yr and ~ 8 .7 x 1 4 yr. Compared to the 
median age of the ONC binaries of 1 Myr, derived from Tab. [6] 
this is very short. These high accretion rates could not have been 
sustained over the entire early evolution process even if the disk 
masses were initially ten times more massive than we now ob- 
serve. If both components were not observed at a younger age 
than assumed or in a short-lived above-average state of accre- 
tion, binary component disks would need substantial replenish- 
ment to display the strong accretion activity we detect. 



4.3. Is the existence of an inner disk linked to planet 
formation in binaries? 

Primordial disks like those detected around the binary compo- 
nents in this sample contain the basic material for the forma- 
tion of planets. Hence, any peculiarities in the evolution of disks 
in these systems can leave their footprints on the population of 
planets in binaries, and it should be instructive to compare the 
properties of planets and disks around the individual components 
of binaries. 

A recent census identified 40 planets in 35 multiple sys- 
tems (Eggenberger & Udry 2010, and references therein), al- 
most all of which orbit the more ma ssive component of 
the binary dMugrauer & Neuhauser 2009). An additional eight 
planets were claimed to re s ide in circumbinary (P-type) or- 



bits (PSR B 1620-2 6. iRasl 



7 Compare also singles in Taurus from Muzerolle et al. ( 1998), which 
are as well at similar accretion rates. 



_ _ J 1994 H P 202206. | Correia et al l 

12005b HW Vir. iLeeet al.ll2009t NN Ser. iBeuermann et al.l l20 10: 
DP Leo. lOian et alJ^OlOafoS Vir. lOian et aLpOlObT HU Aqr, 
lOian et alJl201lh Kepler- 1 6. iDo vie et al.ll201 If) . Although some 
of these planets still need confirmation, we note that all lat- 
ter candidate hosts are spectroscopic binaries with comparably 
small separations. From the thus composed picture of planet oc- 
currence in binaries (~50 planets in multiples, ~2 around the less 
massive component, 8 circumbinary planets), one might suspect 
that (i) planet formation around the less massive components of 
binaries is suppressed and (ii) that circumbinary planets are rare, 
but do exist. These observations might either be caused by selec- 
tion effects, since spectroscopic binaries and fainter secondary 
stars are less often targeted by spectroscopic surveys, or be the 
consequence of peculiar disk evolution in binaries. 

In an attempt to carefully eval uate systematical errors and 
biases, lEggenberger & Udrvl d2010t) discovered an underdensity 
of planets around stars with close, 35-225 AU stellar compan- 
ions when compared to single stars with similar properties. This 
is not seen for wider binaries. The upper limit of 225 AU notice- 
ably coincides with our 200 AU transition to synchronized disk 
evolution. This suggests a common origin of both effects. To test 
whether shorter lived disks in close binaries can explain the defi- 
ciency of planets in binaries of the same separation range, we de- 
rived the disk frequency in both subsamples. Binaries with sepa- 
rations <200 AU have a slightly smaller fraction of disk bearing 
components (34+g%) than wide >200 AU pairs (37+g%). This 
differenc e is, however, not as pronoun ced as the 1. 6-2. lcr differ- 
ence that lEggenberger & Udrvl d2010h observe for the frequency 
of planets in the close and wide sample. We concluded that either 
our sample is not large enough to reveal a significant difference 
or planet formation in binaries is not strongly correlated with the 
occurrence of accretion disks. 

If the apparent paucity of planets orbiting the less massive 
components of binaries were not entirely due to selection effects, 
it could be a consequence of differential disk evolution, which 



18 



Daemgen, Correia, & Petr-Gotzens: Protoplanetary disks of T Tauri binaries in the ONC. 



we observe as mixed systems with only one accreting compo- 
nent. If the probability of a star to eventually host a giant ex- 
trasolar planet were significantly correlated with the lifetime of 
its disk, we would expect mainly CW systems, to evolve into 
circumprimary planetary systems while WC-type binary-disk- 
systems would preferably evolve into circumsecondary systems. 
However, we do not observe any significant difference in the ap- 
pearance of CW systems versus WC (see Sect. 14.2.31 ). In the 
same way, we do not see any binaries in which both compo- 
nents are orbited by their individual planets, although the ma- 
jority of disks evolves s ynchronizedly. In agreem ent with this 
result are the findings of iJensen & Akesonl (120031) . who found 
that disk masses around the primary are always higher than sec- 
ondary disk masses in four T Tauri binaries, independent of their 
classification as CC, CW, WC, or WW. Again, a possible expla- 
nation could be that the evolution of the inner (accretion) disk 
is not strongly related to the formation of planets and that other 
factors dominate the planet formation process. 



5. Summary 

We have presented high-spatial-resolution near-infrared spectro- 
scopic and photometric observations of the individual compo- 
nents of 20 young, low-mass visual binaries in the Orion Nebula 
Cluster. The sample was comple mented with s imilar observa- 
tions of six additional targets from lCorreia et al.ld2012l in prep.). 
We have measured the relative positions, JHK^ photometry, 
and /iT-band spectra including the accretion-indicating Brackett- 
y feature in order to derive the projected binary separations as 
well as the absolute magnitude, spectral type, effective temper- 
ature, extinction, veiling, luminosity, and the probability of dust 
and accretion disks around each binary component. By placing 
the components into an HR diagram and comparing with pre- 
main sequence evolutionary tracks, we have estimated the indi- 
vidual age, mass, and radius, as well as the mass accretion rate 
for each individual component. 

Putting the results into context with the star forming envi- 
ronment of the Orion Nebula Cluster and with other young low- 
mass binary studies, we conclude the following: 

1. We have found evidence of a slightly lower frequency of 
circumstellar disks around the individual components of bi- 
naries compared to that around single stars of the ONC, in 
agreement with theory. We have measured a corrected accre- 
tion disk fraction of 40*i°% for stars in multiple systems of 
the ONC, which is lower than the ~50% a ccretion disk frac- 
tion o f single stars in the ONC found by Hille nbrand et alj 
(1998). A similar result was found for dust disks, as indi- 
cated by NIR excess emission, although with lower signifi- 
cance. As observed for single stars in other clusters, binary 
components of the ONC have been more often found to con- 
tain dust disks than accretion signatures. 

2. The evolution of disks around both components of a binary 
is correlated for binaries with separations of 200 AU and be- 
low. This was inferred from our inability to detect any mixed 
pairs of accreting and non-accreting components with sepa- 
rations <200 AU, and that the populations of mixed pairs ex- 
hibit significantly larger separations (99.5% confidence) than 
pairs of two accreting or two non-accreting components. We 
have demonstrated that this synchronization is probably not 
caused by a feeding mechanism involving a circumbinary 
disk, but possibly instead closer binaries that harbor equal- 
mass components. 



3. Mixed pairs including an accreting primary and those with 
an accreting secondary have been observed to be almost 
equally abundant. In addition to the implication that mixed 
pairs are common, there is apparently no preference for ei- 
ther the disk of the more or less massive binary component to 
dissolve first. This points to a weak correlation between the 
binary mass ratio and the presence of a disk around either 
binary component. 

4. We have found that the mass accretion rates of binary com- 
ponents in the ONC do not differ from the accretion rates of 
single stars and binary components in Orion and Taurus, re- 
spectively. Since disk masses and radii of primordial disks in 
the ONC are - on average - lower, this can potentially lead 
to shorter lifetimes of disks around binary components, in 
agreement with our finding of fewer dust and accretion disks 
than singles of the ONC. 

5. We have measured no strong correlation between the exis- 
tence of planets around the components of main-sequence 
binaries and the occurrence of accretion features mea- 
sured around young binary stars in this paper. Although 
planets seem to be slightly suppressed in binaries of 
separations smaller than ~ 200AU (1.6-2.1cr significance; 
Eggenber ger & Udry||2010h . in wider binaries they are not, 
and we have not found any equally large differences between 
the presence of disks in close and wide binary systems. 
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Appendix A: Statistical evaluation of the disk 
frequency 

We used a Bayesian approach to derive the probability density 
function (PDF) of the disk frequency around the individual com- 
ponents of binaries in the ONC. The location of the maximum 
of the distribution and its width represent the most robust es- 
timates of our disk frequency and its uncertainty, respectively. 
Input information comes from the measured values of the Bry 
equivalent width and the continuum noise measurements in each 
target spectrum. 

A. 1 . Derivation of the formulae 

We defined the probability density P(§\D), given the data D, of 
the true disk frequency § with < # < 1 to be in the interval 
[#,# + d#], as 

P(#\D) = Yprob (#,k=j\D) (A.l) 

7=0 
n 

= 2 prob = j, D) ■ Pn (k = j\D) . (A.2) 

j=o 

While prob(#, k\D) describes the PDF as a function of the disk 
frequency & and the observed number of disks k, the prob(#|fc, D) 
in (IA.2I) are the probability density functions PDFt, which de- 
scribe the probability density of finding the true disk frequency 
to be in the range # + d§] when the number of detected disks 
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is fixed at k (see 9A.1.U . The latter were weighted by p„(k\D) 
in order to account for their relative probability (see jjA.1.21 ). To 
include the contributions of all possible outcomes (0. . . n disks 
surrounding n targets), we summed up all partial PDFs to the 
final PDF P(§\D). 

A.1 .1 . The probability density for k disks in n targets 

The PDF prob(#|A:, D) can be most accurately described as being 
binomial, since there are only two options "disk" or "no disk". 
For instance, starting with Bayes' theorem, we found that 



Table A.l. Accretion 



p(§\k, n) 



p (k\§, n) p (fl|n) 
p(k\n) 

-fff~ k 



1 

17+ 1 



(n + l) 



#) 



n—k 



1 for < ■& < 1 
else 



(A.3) 



(A.4) 



(A.5) 



for the total number of targets n. The prior p(ft\n) was assumed 
to be uniform, since it simply describes the probability of find- 
ing a disk fraction anywhere between and 1 . In addition, we 
found that p(k\n) = l/(n + 1), since it describes the probability 
of selecting a solution among n+l possibilities {0, . . . , n}. 

A.1.2. The weighting factors p(k\D) 

We were able to compare the measurement of the equivalent 
width W(Bry) to the noise limit at the position of the Bry line 
(Table lATTb - Assuming the noise and the underlying distribution 
of the Bry measurement to be approximately Gaussian (which 
can safely be assumed owing to the high photon count rates), we 
were able to derive the probability EL of significant Bry emission 
and hence a disk by integrating a Gaussian distribution (with 
median p. = W(Bry) and width <x = AW(Bry)). The lower and 
upper integration limits were the continuum noise and infinity, 
respectively. For example, if we measured an equivalent width 
of 1.5+0.5 A and the continuum noise level is 1.0 A (i.e. right at 
the lcr limit of our measurement), then the derived probability 
was 11 = 0.84, i.e. we assigned a 16% (=lcr) probability that the 
system contained no disk. 

With this knowledge, we were able to derive p„(k\D), i.e. the 
probabilities of finding k disks in a set of n targets. We defined 
pairs {1 — II,-, n, } of the probabilities of finding "no disk" or "one 
disk" around a particular target i. The p n (k\D) were then derived 
by convolution (zero padded at the edges) of all (1 - fl,-, II,} 



p„(fc|D) e {1 - IIo, IIo} *{!-!!!, I!, 



{i-n„,n„} .(A.6) 



For illustrative purposes, the probabilites of finding 0, 1, or n 
disks in a sample of n was 

n 

Pn (0\D) = (i - no) (i - no ... (i - n„) = \~[ (1 - n,-) (A.7) 

1=0 



Name 


Component 


[A] 


Wmin 

[A] 


[AD95] 1468 


A 


1.48 ±0.49 


1.47 




B 


0.68 ±0.56 


1.62 


[AD95] 2380 


A 


1.00 ±0.19 


0.65 




B 


0.85 ± 1.12 


3.51 


JW235 


A 


5.13 ±0.47 


1.21 




B 


2. 14 ±0.54 


1.59 


JW260 


A 


-2.11 ±0.13 


0.37 




B 


-4.05 ±0.15 


0.48 


JW553 


A 


0.02 ±0.14 


0.44 




B 


-0.26 ± 0.24 


0.75 


JW566 


A 


0.63 ±0.1 8 


0.53 




B 


0.45 ±0.23 


0.70 


JW598 


A 


-0.33 ±0.10 


0.34 




B 


1.10 ±0.66 


2.06 


JW648 


A 


0.42 ±0.17 


0.49 




B 


0.48 ±0.31 


0.84 


JW681 


A 


-0.38 ±0.13 


0.46 




B 


0.94 ±0.21 


0.68 


JW687 


A 


0.09 ±0.1 6 


0.43 




B 


-0.05 ±0.16 


0.51 


JW876 


A 


0.52±0.18 


0.54 




B 


0.96 ±0.18 


0.51 


JW959 


A 


-1.48 ±0.22 


0.65 




B 


-0.82 ±0.25 


0.76 


[HC2000] 73 


A 


1.40 ±0.14 


0.41 




B 


-0.10 ±0.49 


1.35 


TCC15 


A 


-0.79 ±0.21 


0.67 




B 


16.69 ±2.73 


5.44 


TCC52 


A 


1.51 ±0.10 


0.28 




B 


1.14 + 0.14 


0.40 


TCC55 


A 


-0.14 ±0.20 


0.63 




B 


-0.24 ± 0.23 


0.68 


JW63 


A 


-0.04 ±0.34 


1.02 




B 


0.61 ±0.43 


1.17 


JW128 


A 


0.53 ±0.49 


1.28 




B 


0.64 ± 0.49 


1.28 


JW176 


A 


0.34 ± 0.47 


1.26 




B 


0.22 ±0.58 


1.47 


JW391 


A 


7.87 ±0.26 


0.69 




B 


1.41 ±0.48 


1.32 


JW709 


A 


-0.16±0.36 


1.01 




B 


0.42 ±0.35 


0.97 


JW867 


A 


0.91 ±0.27 


0.69 




B 


2.37 ±0.29 


0.77 



Notes. (a) The equivalent width as measured in our spectra. No cor- 
rection for veiling was applied, since W m j n and Wbi T scale the same 
with veiling. {h] Equivalent width equivalent to 3x continuum noise 
level around Bry. This is the assumed detection limit for significant 
Bry emission. 



Pn (i\D) = 2in (i_n ^ ) n ' ' 

1=0 J 
Pn (n\D) = Y\nj . 

7=0 

Other p n (k\D) could take more complicated forms. 



(A.8) A.2. Application to the ONC 

For all calculations, we assumed that n = 42, which is equal 
(A 9) to tne numrjer °f target components excluding JW 260 (see 
Sect. 14.2. U . The individual n,- were calculated from our Bry 
measurements in Tab. IA.ll The final probability density was de- 
rived from ( 1A.2I ) and is depicted in Fig. IA.ll 
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Fraction of disks in binary components 



Fig.A.l. The resulting ONC probability density function as a 
sum of the individually weighted p(ft\k, D). The resulting best- 
fit value of the disk fraction of binaries in the ONC (see bar) is 
p _ o ^s +0 09 

r — u.JJ_ 00g . 

We found a best fit value of F = 35*g % for the disk frequency 
in the ONC binaries. The uncertainties of F are defined by the 
16% (= lcr) wings of the distribution. 
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Table 6. Individual component properties of all targets with detectable photospheric features 









7dr° 


Av 




L, 


age' J 


M," 


R, b 


Ej-tf 






disk 






Name 


Comp. 


SpT 


IK] 


[mag] 


(F Ke JF K ,) 




[Myr] 


[M s ] 


[Re] 


[mag] 


[mag] 


[A] 


prob. 1 ' 


log(L acc /L ) 


[lO-'Moyr 1 ] 


[AD95] 1468 


A 


M4± 1 


3270 ± 150 


3.7 ±1.1 


0.00 ± 0.29 


0.11 ±0.03 


4.0 ±« 

1.4 


0-23 ±™ 

0.05 


1.04 ±0.05 


0.43 


0.47 


1.48 ±0.50 


0.50 


-1.92±0.21 


2.2 ± 1.44 




B 






0.5 ±1.3 














(<1.62) 


0.05 






[AD95] 2380 


A 
B 






















(1.00±0.19) 
(<3.51) 


0.97 
0.01 






JW 235 


A 






0.0 ± 1.0 
















(5.13 ±0.47) 


1.00 








B 


M3±3 


3420 ±440 


5.7 ±0.9 


0.73 ±0.65 


0.62 ±0.16 


i.o ±11 


35 ± 23 


2.24 ± 0.29 


-0.04 


-0.26 


3.67 ±0.93 


0.84 


-1.37 ±0. 18 


11 ±g 


JW 260 


A 


G0± 1 


6030 ± 170 


0.0± 1.1 


0.20 ± 0.04 


18.0±5.5 


3.5 ±0.5 


2 5± a2 


3.89 ±0.11 


0.29 


0.19 


-2.52±0.15 


0.00 








B 


F7±3 


6320 ±210 


0.0 ±1.1 


0.00 ±0.04 


12.1 ±3.7 


4.0 ±0.5 


2 3 ± °- 2 
i - J 1 0.2 


2.91 ±0.10 


0.44 


0.43 


-4.05 ±0.15 


0.00 






JW553 


A 


K7± ! 


4060 ±250 


3.0 ±0.1 


0.22 ±0.06 


5.2± 1.5 


0.35±Sjl 


0.8 ±i" 


4.63 ±0.28 


-0.01 


-0.03 


<0.54 


0.00 


< -1.11 


< 18 




B 


Ml ±3 


37 10 ±440 


1.4 ±1.7 


1.72 ±0.57 


0.16 ±0.02 


8±F 


0.57 ±« 3 


0.97 ±0.12 


-0.14 


0.26 


<2.04 


0.00 


< -2.68 


<0.14 


JW566 


A 


K7± 1.5 


4060 ± 280 


4.4 ± 1.1 


0.02 ± 0.05 


1.3±0.2 


1 2± LS 

5 


0.75 ±9* 


2.33±0.16 


0.49 


0.78 


0.64 ±0.19 


0.70 


-1.01 ±0.05 


12 ± 5 




B 


M1.5±2 


3705 ±290 


0.0 ±0.3 


0.00 ±0.05 


1.2±0.1 


9 ± 0:5 


48 ± 2 


2.61 ±0.20 


0.21 


0.29 


<0.70 


0.13 


< -1.30 


< 11 


JW 598 f 


A 
B 


K5±2 
M2±4 


4350 ±380 
3560 ±570 


4.4 ± 0.4 
5.0 ±0.4 


0.01 ±0.44 










(<0.34) 
<2.08 


0.00 
0.07 






JW648 


A 


M0± 1 


3850 ± 170 


2.3 ±0.1 


0.06 ±0.05 


1.7 ±0.1 


0.7 ±0.2 


57 ± 013 
0.30 ±«J 


2.90 ±0.13 


0.16 


0.30 


<0.52 


0.33 


< -1.43 


<7.6 




B 


M3.5±2 


3340 ± 290 


2.3 ±0.3 


0.00 ±0.14 


0.49 ± 0.04 


1.0 ±0.8 


2.09 ±0.18 


0.21 


0.14 


<0.84 


0.12 


< -1.76 


< 4.9 


JW681 


A 
B 


Ml ± 1.5 
M3.5±4 


3710±220 
3340 ± 550 


8. 1 ± 0.6' 
7.6 ±1.2/ 


0.09 ±0.06 
0.26 ±0.15 












<0.50 
1.18 ±0.26 


0.00 
0.90 






JW687 


A 


M2.5±°- 5 


3490 ± 150 


6.1 ±0.1 


0.43 ±0.10 


1.7±0.1 


0.10 ±0.05 


0.38 ±™ 


3.58 ±0.15 


0.03 


-0.06 


<0.61 


0.02 


< -1.56 


< 10 




B 


M2± 1 


3560 ±150 


7.4 ±0.4 


1.62 ±0.27 


0.54 ±0.06 


1 3 ± 05 

0.4 


39 ± 05 

U x 0.06 


1.93 ±0.08 


0.34 


0.49 


<1.34 


0.00 


< -1.67 


< 4.2 


JW 876 


A 
B 


M0.5 ± 1 
Ml. 5 ± 3 


3790 ±210 
3630 ±440 


2.3 ±0.6' 
0.0 ± 1.0^ 


0.02 ± 0.06 
0.38±0.13 










<0.55 
1.30±0.24 


0.47 
0.99 


< -0.87 
-0.97 ± 0.04 




JW 959 


A 
B 


K3±3 
K3 ± 2 


4730 ±450 
4730 ±370 


3.0 ±0.8/ 
5.6 ±0.7^ 


0.27 ±0.10 
0.00 ± 0.08 














-1.86 ±0.27 
-0.81 ±0.25 


0.00 
0.00 






[HC2000] 73 


A 


M2± 1 


3560 ±150 


1.1 ±0.5 


1.63 ±0.28 


0.21 ±0.03 


3.5 ±|;° 


038 *m 


1.21 ±0.05 


0.34 


0.51 


3.68 ±0.37 


1.00 


-1.78±0.16 


2 1 ± 1 1 




B 


M7± 1 


2880 ±140 


0.0 ± 0.5 


0.01 ±0.24 


0.08 ±0.01 


0.6 ±« 


09 ± 

0.05 


1.17 ±0.06 


-0.22 


-0.02 


<1.36 


0.00 


< -2.64 


< 1.2 


TCC 15" 


A 


K4±2 


4590 ±290 


10.6 ±0.4 




2.7 ±0.3 


1.8 ±« 


1.5 ±°; 2 


2.63 ±0.17 


0.09 


0.06 


(-0.79 ±0.21) 


0.00 








B 






3.0 ±2.8 












(16.69 ±2.73) 


1.00 






TCC52 


A 


M0.5 ± 2 


3790 ± 300 


0.7 ±0.3 


4.42 ±0.36 


10.4 ±0.9 


01 ± 020 


0.6 ±g| 


7.50 ±0.59 


0.36 


0.59 


8.00 ±0.56 


1.00 


0.39 ±0.20 


1230 ±|° 




B 


M2± 1 


3560 ±150 


2.0 ±0.4 


0.80±0.12 


3.1 ±0.4 


0.05 ±111 


39 ± 

U- ^ 1 0.05 


4.61 ±0.19 


0.28 


0.40 


2.04 ±0.25 


1.00 


-0.49 ±0.07 


153 ±« 


TCC 55 


A 


M3± 1 


3420 ±150 


3.0 ±0.7 


0.14 ±0.10 


0.09 ± 0.02 


6±" 2 


29 ± °- 07 

X 0.06 


0.87 ±0.04 


0.86 


1.23 


<0.72 


0.00 


< -2.19 


< 0.7 




B 






0.0 ±0.7 














(<0.68) 


0.00 






JW63 


A 


K7±2 


4060 ± 320 


2.0 ±0.3 


0.00 ±0.05 


1.3 ± 0.1 


1 mZ - 0.5 


u./j 2: 25 


2.32 ±0.18 


-0.06 


-0.10 


<1.02 


0.00 


< -1.54 


< 3.6 




B 


M3 ± 2 


3420 ± 290 


3.2 ±0.3 


0.00 ±0.21 


0.84 ±0.07 


Q7 + 
u - ' 1 0.6 


r\ •)•) , f) 1 5 

°- 33± a!o 


2.62 ±0.22 


-0.16 


-0.44 


<1.17 


0.09 


< -1.76 


< 5.5 


JW 128 


A 


MO ±2 


3850 ±390 


2.2 ±0.2 


0.00 ±0.11 


1.6 ±0.1 


7 ± 1 

, *B 


0.56 ±°|| 


2.81 ±0.28 


-0.05 


-0.05 


<1.28 


0.06 


< -1.21 


< 12 




B 


MO ±2 


3850 ±390 


0.8 ±0.2 


0.00 ±0.1 3 


0.89 ±0.06 


1.2 ± 2 j 


r\ c-i 40 
x 0.19 


2.12 ± 0.21 


0.07 


0.14 


<1.26 


0.09 


< -1.40 


< 5.9 


JW 176 


A 


K7±2 


4060 ± 320 


3.2 ±0.2 


0.02 ±0.16 


3.5 ±0.2 


45 ± 1145 

W ' HJ x 0.40 


0.75 ± 


3.80 ±0.30 


-0.04 


-0.08 


<1.50 


0.02 


< -0.81 


< 31 




B 


MO ±2 


3850 ±390 


4.0 ±0.2 


0.00 ±0.1 9 


2.8 ±0.2 


0.45 ± 0.40 


58 ± ™ 


3.75 ±0.38 


-0.18 


-0.26 


<1.47 


0.02 


< -0.93 


< 30 


JW391 


A 


M0.5 ± 3 


3790 ±440 


1.2±0.2 


0.71 ±0.31 


1.1 ±0.1 


9± 1 6 


52 ± 


2.42 ±0.28 


0.21 


0.35 


12.91 ±0.45 


1.00 


-0.12±0.11 


140 ± 120 




B 


M3.5±2 


3340 ±290 


0. 1 ± 0.2 


0.00 ±0.22 


0.15±0.01 


3 5 ± 35 
JJ x 1.5 


0.28 ±g| 2 


1.16 ±0.10 


0.20 


0.11 


1.41 ±0.48 


0.58 


-2.20 ±0.22 


1 1 ± °- 8 

1 ■ 1 ^ 0.8 


JW 709 


A 


K7±2 


4060 ±320 


1.7 ±0.2 


0.00 ±0.24 


0.94 ±0.06 


1.8 ±" 


0.76 ± 040 


1.97 ±0.16 


-0.00 


-0.01 


<1.01 


0.00 


< -1.66 


< 2.3 




B 


K7+ 1 


4060 ±250 


1.6±0.2 


0.00 ±0.11 


0.69 ±0.05 


3.0 ±« 


0.78 ±|| 


1.68±0.10 


-0.01 


-0.02 


<0.97 


0.06 


< -1.87 


< 1.2 


JW 867 


A 


Ml ± 1 


3710 ± 150 


0.1 ±0.3 


0.13 ±0.12 


0.89 ±0.08 


1.0 ± 2- 4 
9 ± 


048 *m 


2.29 ±0.09 


0.16 


0.24 


1.03 ±0.31 


0.80 


-1.50±0.12 


6.0 ±2.3 




B 


M2.5 ± 1 


3490 ±150 


0.0 ±0.3 


0.02 ±0.11 


0.77 ±0.07 


0-37 *gg 


2.41 ±0.10 


0.14 


0.19 


2.41 ±0.29 


1.00 


-1.05 ±0.05 


23 ±4.8 



Notes. <D) Effective temperatures were conve rted with the temperat ure scale for pre-main sequence stars given in iLuhman et alj |2003) for spectral types later than M0. The effective 
temperatur e s of e arlier spectral types are from lSchmidt-KaleJ (l982). Uncertainties are from r c(f at the SpT uncertainty limits. (A) Results from comparing T sff and L, with the models of 
Siess et al. (2000). w Positive values indicate emission lines, negative values indicate Bry in absorption. Upper limits refer to the minimum equivalent widths of an emission line at A Bry 
that can be detected with 3<r significance. Whenever possible, equivalent widths were extracted from the r^-corrected spectra. Values from targets without r^are in parantheses. These can 
be regarded as lower limits to the equivalent width of the Bry line and no accretion luminosities or mass accretion rates were derived. The disk probability is the result of the integration 
described in 9A.1.2I This does not exclude the possibility of 'filled-up' absorption lines where emission and absorption cancel each other. (c) The spectral features allow us to assign a 
spectral type, but the spectral fitting routine was unable to match the overall continuum shape of this target, in particular for a strong increase at A > 24000 A. These extinction values 
were derived from our spectral fitting instead of photometry, since our photometric information was incomplete. 



